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ABSTRACT 


A  microprocessor-based  digital  flight  control  system  (Micro- 
DFCS)  is  implemented  in  Princeton's  Variable-Response  Research  Air¬ 
craft  (VRA),  and  evaluated  through  flight  testing.  The  flight  control 
computer  program,  CAS-1,  provides  three  longitudinal  control  options; 
direct  (unaugmented)  command,  pitch  rate  command,  and  normal  accelera¬ 
tion  command.  The  latter  two  options  are  Type  0  systems  designed  by 
linear-quadratic  control  theory.  The  flight  tests  are  designed  to 
investigate  the  characteristics  of  these  control  laws  in  flight,  to 
provide  experimental  evidence  for  sampling  rate  requirements,  and  to 

demonstrate  digital  flight  control  of  the  VRA  with  the  Micro-DFCS. 

* 

The  Micro-DFCS  is  installed  in  the  VRA  and  functions  properly. 
The  qualitative  and  quantitative  flight  tests  indicate  that  sampled- 
data  regulator  theory  adequately  accounts  for  the  phase  lags  associated 
with  sampling  and  holding  continuous  signals.  Low  sampling  rates,  of 
10  samples  per  sec  (sps)  or  less,  can  be  used  as  a  consequence.  The 
controllability  and  responsiveness  resulting  at  each  of  the  three 
command  modes  are  good  at  10  sps.  The  pitch  rate  command  option  is 
considered  to  give  good  response  down  to  7  sps,  although  the  vibrations 
generated  by  the  discrete  elevator  movements  sound  somewhat  like  the 
stall  buffet,  which  is  disconcerting  to  the  pilot.  For  all  sampling 
rates  tested  (except  5  sps),  pitch  rate  and  normal  acceleration  control 
improves  response  over  the  direct  command  option;  however,  the  discrete 


movements  of  the  elevator  are  more  apparent  for  7  sps  and  below. 

Experimental  results  indicate  that  the  Micro-DFCS  provides 
a  practical  way  to  implement  digital  flight  control.  More  research 
is  necessary  to  investigate  digital  flight  control  using  advanced 
command  and  control  modes. 
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INTRODUCTION 


A  properly  implemented  Command  Augmentation  System  (CAS) 
provides  stability  and  control  augmentation  which  can  significantly 
improve  the  basic  aircraft  handling  qualities.  A  digitally  mechanized 
CAS,  as  opposed  to  one  which  uses  analog  components,  allows  advanced 
concepts  of  control  theory  to  be  easily  reduced  to  practice.  Using 
the  computational  capabilities  of  digital  computers,  advanced  control 
modes  can  be  implemented  and  scheduled  according  to  the  task  and  flight 
condition. 

Digital  microprocessors  and  modern  control  theory  are  two 
relatively  new  areas  for  research  in  the  Digital  Flight  Control  (DFC) 
context.  The  application  of  both  technologies  to  DFC  will  help  realize 
the  full  potential  of  flight  control.  To  date,  few  aircraft  have  been 
flown  under  digital  control  and  none  have  used  a  microprocessor  based 
digital  flight  control  system  (Micro-DFCS)  for  a  manned  vehicle.  Like¬ 
wise,  the  use  of  modern  control  theory  for  DFC  in  actual  flight  tests 
has  been  limited.  More  flight  testing  of  DFC  using  these  concepts 
is  necessary  to  fully  understand  the  relationships  between  theory  and 
practice. 

This  thesis  describes  an  experimental  program  for  investigating 
DFC  using  a  Micro-DFCS.  Simple  Type  0  longitudinal  control  laws,  de¬ 
signed  with  modem  control  theory,  are  implemented  on  a  Micro-DFCS 
installed  in  the  research  aircraft.  Actual  flight  tests  are  conducted 
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to  provide  results  and  experimental  evidence  for  DFC  and  modern  control 
theory. 

1.1  ADVANTAGES  OF  DIGITAL  FLIGHT  CONTROL  USING  MODERN  CONTROL  THEORY 

Modern  control  theory  combines  state  space,  time  domain,  and 
optimal  control  concepts  with  earlier  frequency  domain  methods  in  the 
design  of  control  systems.  The  procedure  uses  high-speed  digital 
computers  to  generate  control  gains  and  parameters.  This  method  allows 
easy  development  of  advanced  control  modes  such  as  controllers  with 
state  and  control  rate  restraint,  control  decoupling,  interconnects, 
and  Type  1  CAS  structures  (although  this  study  examines  none  of  these 
advanced  modes).  For  a  Digital  Flight  Control  System  (DFCS),  it  is 
important  to  minimize  sampling  rates  to  reduce  computational  require¬ 
ments.  High  sampling  rates  will  be  required  unless  care  is  taken  in 
the  modeling  and  design  process  (Ref.  1);  hence,  modem  digital  control 
design  techniques,  including  both  time-domain  and  frequency-domain 
methods,  generally  model  the  sampling  effects  well  enough  to  allow  low 
rates  to  be  used  without  loss  of  stability  or  performance. 

Some  of  the  major  advantages  of  using  a  digital  computer 
for  flight  control  are  the  following: 

•  Ease  of  scheduling  gains  and  changing  control  modes 

•  Ability  to  implement  advanced  control  modes  easily 
(Ref.  2) 

•  Ease  of  changing  performance  characteristics 

•  Ability  to  time  share  logic  to  serve  a  number  of  other 
functions  (such  as  navigation  or  fuel  management)  (Ref. 3) 

•  Ability  to  perform  in-line  self  testing,  resulting  in 
improved  performance  (Ref.  2) 
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The  first  advantage  listed  is  attributed  to  the  digital  computer's 
conditional  logic.  By  monitoring  the  aircraft's  flight  condition  and 
pilot's  requests,  the  computer  can  make  decisions  to  change  the  param¬ 
eters  of  the  presently  running  control  law  or  can  implement  a  completely 
different  CAS  structure.  The  last  four  points  listed  are  due  mainly  to 
the  flexibility  of  designing  with  software.  Once  the  pilot's  commands 
and  the  feedback  variables  are  interfaced  with  the  computer's  analog 
input  channels,  the  software  can  be  developed  to  manipulate  the  input 
data  in  a  variety  of  ways  to  perform  the  desired  functions.  Advanced 
control  modes  would  require  different  sets  of  calculations  i:i  software, 
but  no  increase  in  hardware  complexity.  Performance  characteristics 
of  the  DFCS  can  be  changed  by  simple  software  modifications  instead 
of  a  complete  redesign  of  hardware.  The  computer  can  be  programmed  to 
diagnose  malfunctions  in  the  control  system  and  take  corrective  action. 

For  instance,  the  computer  can  compare  the  output  from  feedback  sensors 
with  the  values  expected,  given  by  a  state  estimator.  If  a  sensor's 
value  is  completely  different  from  the  expected  value,  a  back-up  sensor 
can  be  brought  on  line. 

It  is  becoming  increasingly  apparent  that  the  principal  means 
of  implementing  the  DFCS  in  the  future  will  be  with  microprocessor 
technology.  Great  advances  within  the  last  few  years  in  the  development 
of  large-scale-integrated  circuits  enable  these  small,  inexpensive 
components  to  perform  the  arithmetic  logic,  data  storage,  control  se¬ 
quencing,  and  input/output  functions  necessary  for  the  DFCS.  A  Micro-DFCS 


will  give  a  significant  reduction  in  cost,  size,  weight,  and  power 
usage  for  the  flight  control  hardware.  These  four  characteristics  will 
probably  have  the  most  direct  impact  on  cost /performance  criteria  and 
therefore  will  be  the  immediate  stimulus  for  using  the  Micro-DFCS. 

The  present  major  limitation  of  microprocessors  is  the  small  data  word 
length.  Most  microprocessors  on  the  market  today,  have  word  lengths 
of  four  or  eight  bits  (although  a  few  have  16  bit  words).  This  re¬ 
duces  the  computational  speed  relative  to  minicomputer  systems,  increas¬ 
ing  the  need  for  low  sampling  rates. 

Very  few  flight  test  programs  have  been  done  to  research 
digital  flight  control.  To  date,  the  tests  used  aircraft  and  experi¬ 
mental  equipment  with  high  operating  costs  (due  to  their  large  size, 
high  fuel  use,  and  complexity).  References  2  and  4  to  10  describe 
flight  test  programs  that  used  a  DFCS.  Because  of  the  high  cost  of 
these  programs,  the  data  on  actual  implementation  of  digital  flight 
control  is  understandably  limited.  No  flight  tests  prior  to  this 
study  have  used  microprocessor  technology  for  digital  flight  control 
in  manned  vehicles. 

Similarly,  the  availability  of  flight  test  results  for  a  DFCS 
designed  with  modern  control  theory  is  small  (References  4,  S_,  10  and 
11  document  investigations  that  have  used  modem  control  laws  in  flight, 
and  there  may  be  other  examples  not  cited  here).  Past  programs  have 
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usually  used  the  conservative  approach  of  digitizing  an  analog  control 
law  which  leads  to  high  sampling  rates  (usually  50  per  sec  or  higher). 
Sampled-data  regulator  theory  (Ref.  12)  gives  control  law  designs  that 
take  into  account  the  discrete-time  nature  of  the  calculations  and 
makes  more  efficient  use  of  digital  systems  by  allowing  lower  sampling 
rates.  More  experimental  data  is  required  through  flight  testing  to 
define  the  methodologies  and  verify  the  concepts  of  digital  flight 
control  using  microprocessor  technology  and  modern  control  theory. 

Such  testing  will  promote  the  practical  acceptability  of  digital  flight 
control  and  modern  control  theory  in  the  future. 

1.2  RESEARCH  OBJECTIVES  AND  SUMMARY  OF  MAJOR  RESULTS 

The  flight  testing  of  this  investigation  is  conducted  using 
Princeton's  Variable-Response  Research  Aircraft  (VRA)  operating  from 
the  Flight  Research  Laboratory  (FRL) .  See  Appendix  C  for  a  description 
of  the  VRA  and  FRL.  The  analog  Fly-by-Wire  (FBW)  system  of  the  VRA  is 
augmented  with  the  Micro-DFCS  described  in  Chapter  3  and  Appendices  D 
and  E.  Qualitative  and  quantitative  data  are  collected  on  the  perform¬ 
ance  of  the  system  operating  with  simple  Type  0  control  laws  where  the 
optimal  gains  are  not  scheduled.  Pitch  rate  and  normal  acceleration 
CAS  structures  with  two-state  feedback,  one  control  input  to  elevator, 
and  one  pilot  input  of  longitudinal  stick  are  tested. 

The  testing  is  constructed  to  allow  three  primary  objectives 
to  be  met : 
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•  Provide  results  for  the  investigation  of  optimal  control 
laws  in  flight 

•  Provide  experimental  evidence  for  sampling  rate  require¬ 
ments 

•  Successfully  demonstrate  the  VRA  FBW  system  augmented 
with  the  Micro-DFCS 

The  first  and  second  objectives  are  aimed  at  answering  some  basic  questions 
about  implementing  a  Micro-DFCS.  They  are  accomplished  by  setting  a  varia¬ 
tion  of  sampling  rates  and  optimal  gains  for  the  pitch  rate  control 
mode,  the  normal  acceleration  control  mode,  and  a  direct  (open-loop) 
mode.  The  different  control  configurations  are  evaluated  in  the  landing 
approach,  in  a  tracking  task  at  altitude,  and  by  examining  the  step  response 
of  the  commanded  state  to  a  step  input  with  longitudinal  stick.  The 
telemetry  records  and  pilot's  comments  for  each  variation  are  examined 
to  demonstrate  flying  qualities  and  to  define  advantages  and  limitations 
of  the  Micro-DFCS  using  optimal  control  laws.  The  test  results  are  used 
also  to  establish  some  lower  limits  on  sampling  rate  for  these  particular 
control  configurations.  The  third  objective  is  intended  to  demonstrate 
that  the  Micro-DFCS  equipment  and  related  software  function  in  actual 
flight  as  designed.  Of  course,  the  successful  completion  of  the  first 
two  objectives  depend  on  accomplishment  of  the  third  objective. 

Many  results  were  obtained  from  the  development  and  testing  of 
the  Micro-DFCS;  a  few  major  ones  are  summarized  here.  The  hardware  of 
the  Micro-DFCS  is  completely  installed  and  working.  No  hardware  problems 
were  encountered  except  for  a  structural  vibration  that  needed  to  be 
filtered  out  of  the  normal  accelerometer.  A  sampling  rate  of  10  samples  per 


second  (sps)  is  very  acceptable  for  the  pitch  rate  control  laws  studied. 

At  sampling  rates  of  7  sps  and  below,  the  step  movements  of  the  elevator 
are  very  annoying  and  sound  somewhat  like  a  stall  buffet.  The  present 
results  indicate  that  for  any  given  sampling  rate,  the  pitch  rate  CAS 
improves  the  response  over  that  obtained  with  open-loop  control.  Sam¬ 
pling  rates  of  3  sps  and  above  are  acceptable  for  the  flare  maneuver 
of  the  landing,  but  the  short  final  approach  segment  requires  10  sps 
for  good  control.  Telemetry  records  show  that,  as  expected,  decreasing 
sampling  rates  do  not  affect  the  step  response  of  the  commanded  state 
when  the  gains  are  calculated  using  sampled-data  regulator  theory. 

1,3  ORGANIZATION  OF  THE  THESIS 

The  thesis  is  organized  into  five  chapters  and  five  appendices. 

In  Chapter  2,  the  control  laws  are  developed  with  linear-optimal  control 
theory.  The  optimal  gains  are  calculated  by  a  linear-quadratic  weighting 
method  and  an  implicit  model-following  technique.  The  response  of  the  aug¬ 
mented  system,  as  predicted  by  the  analytical  study,  is  presented  for 
different  control  laws,  sampling  rates,  and  gains.  In  Chapter  3,  the 
basic  aspects  of  the  Micro-DFCS  equipment  and  software  are  outlined  to 
show  how  the  control  laws  are  actually  implemented.  Results  from 
hybrid  tests  (where  the  VRA  is  simulated  on  an  analog  computer) 
are  presented  and  compared  to  the  analytic  results  of  the  previous 
chapter. 

Chapter  4  gives  results  and  analysis  of  the  actual  flight  tests. 


The  theory  and  designs  discussed  in  Chapters  2  and  3  are  evaluated  in 
the  actual  environment  for  which  they  are  intended.  Comparisons  with 
analytic  and  hybrid  test  results  are  made.  Chapter  5  draws  conclusions 
from  this  work  and  makes  recommendations  for  further  study.  Appendices 
A  through  E  provide  details  of  the  VRA  and  FRL  facilities,  a  derivation 
‘of  digital  weighting  matrices,  instructions  for  using  the  APL  functions 
for  generating  optimal  gains,  and  a  detailed  description  of  the  Micro- 
DFCS  equipment  and  related  software. 


2. 


DEVELOPMENT  OF  THE  CONTROL  LAW 


This  chapter  describes  the  optimal  control  theory,  control 
design  objectives,  and  analytical  results  for  the  Micro-DFCS.  Pitch 
rate  and  normal  acceleration  command  control  laws  are  developed  with 
a  second-order  reduced  model  of  the  VRA.  An  easy  extension  of  the 
linear- quadratic  (LQ)  regulator  methods  used  here,  to  dynamic  model 
matching  is  presented. 

2.1  DERIVATION  OF  THE  TYPE  0  DIGITAL  CAS 

The  development  of  the  digital  CAS  in  this  report  uses  sampled- 
data  regulator  theory.  In  this  case  the  selected  control  law  drives  a 
continuous  system  using  piecewise  -  constant  inputs  that  change  only 
at  the  sampling  intervals  (zero-order  hold).  The  objective  of  the 
sampled-data  regulator  is  to  control  the  system  to  be  as  close  as 
possible  to  the  trajectory  obtained  by  using  a  continuous  controller 
(Ref.  12).  The  derivation  proceeds  as  follows. 

The  design  of  the  digital  CAS  starts  with  specification  of 
the  continuous-time  system  dynamics.  Neglecting  disturbance  inputs, 
the  non-linear  equations  of  motion  for  the  flight  vehicle  (Ref.  13) 
can  be  written  as  a  vector  differential  equation: 

x(t)  =  f[x(t),  u(t)]  (1) 
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The  state  vector  x.(t)  is, in  general, of  length  n  and  specifically  contains 
three  components  each  of  translational  rate,  angular  rate,  and  attitude. 


T 

x  *  [u  v  w  p  q  r  6  H]  (2) 

where  (u,  v,  w)  are  body-axis  velocities,  (p,  q,  r)  are  body-axis  rota¬ 
tional  rates,  and  (  0,  <p,  tp  )  are  Euler  angles  of  the  body-axes  with 
respect  to  an  Earth-fixed  frame.  The  control  vector  u(t)  is  of  length 
m  and  includes  deflections  of  the  elevator,  throttle,  flaps,  ailerons, 
rudder,  and  side  force  panels: 

ii  =  I5E  5T  6F  5R  6SF  J 

In  order  to  linearize  the  equations  for  use  in  designing  the  control 
laws,  it  is  necessary  to  divide  the  states  and  controls  into  nominal 
and  perturbation  components, 

*(t)  =  ^(t)  *  Ax(t)  (4) 

H.(t)  =  ^(t)  +  Au(t)  (5) 


The  linearization  is  accomplished  by  taking  a  first-order  Taylor  series 
expansion  of  Eq.  1: 


(6) 


where 


*o(t)  =  l  teo(t)'  J*o[t)  J 

Ax(t)  =  F(t)  Ax(t)  +  G(t)  Au(t) 


(7) 
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The  F  matrix  is  (n  x  n)  and  will  be  referred  to  as  the  system  dynamics 
matrix.  G  is  (n  x  m)  and  is  called  the  control  effectiveness  matrix.  The 
calculation  of  is  still  a  non-linear  problem,  so  the  CAS  will  be 
designed  with  one  specific  x^  in  mind.  F(t)  and  G(t)  are  functions  of 
and  uq;  as  long  as  the  vehicle  flies  within  close  proximity  to  the 
flight  condition  specified,  F(t)  and  G(t)  change  very  little  (Ref.  13). 
This  leads  to  control  laws  based  on  linear  time-invariant  models  which 
could  be  adapted  to  changing  flight  conditions  by  gain  scheduling. 

The  linear  differential  equations  of  motion  (Eq.  7)  are  con¬ 
verted  to  difference  equations  which  take  the  form. 


A4+i  = 


(10) 


where  $  is  the  aircraft's  state  transition  matrix  and  T  is  the  discrete 
control  effectiveness  matrix  (Ref.  12), 


4>  * 


FT 


I  +  FT  + 


(11) 


r  S  F_1 (I  -  $_1)  G 


(12) 


where  T  is  the  sampling  interval. 

To  complete  the  state  space  representation  of  the  aircraft's  dynamics, 
the  output  of  the  system  is  defined  as: 


A*k  =  Hx  %  *  Hu  A“k 


(13) 
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The  Type  0  control  law  takes  the  form  , 


Au^  -  =  -  C  (Ax^-A£  )  (14) 

where  a  starred  vector  represents  the  steady-state  value  of  that  vector 

when  the  input  to  the  CAS  commanded  by  the  pilot  is  attained.  Examining 

Eq.  14,  it  is  seen  that  when  Ax^,  reaches  its  steady-state  value,  the 

output  of  the  controller,  Au^.  ,  will  be  equal  to  the  proper  steady- 

*  * 

state  deflections.  The  steady-state  vectors  Ax  and  Au  are  functions 
only  of  the  pilot's  commanded  inputs  and  are  independent  of  what  feed¬ 
back  gains  are  used  or  how  the  loop  is  closed.  They  are  simply  the 
calculated  values  of  Ax^  and  Au^,  for  a  particular  commanded  input. 

The  states  or  controls  to  be  commanded  by  the  pilot  are  speci¬ 
fied  in  the  control  law  by  defining  values  for  Hx  and  Hu  of  the  output 
equation,  Eq.  13  (Ref.  14).  In  this  way,  the  pilot's  inputs  are  defined 
as  the  system's  outputs  to  be  controlled  by  the  CAS.  In  other  words,  the 
aircraft's  motions  are  the  variables  commanded  by  the  pilot's  controls 
rather  than  the  control  surface  positions.  Since  the  pilot's  inputs 
will  be  considered,  the  system's  outputs  after  the  selection  of  and 
Hu,  Eq.  13  can  be  rewritten , 


A*d, 


=  Hx  A*k  +  Hu  AHk 


(15) 


where  A^ 


contains  the  pilot's  commanded  inputs. 


The  relations  between  A^ 


* 

Ax  ,  and  Au 


* 


are  described  by 


Eqs.  7  and  15  when  Lx  -  0 


defining  the  trim  condition: 


f-  - 
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(16) 


If  the  number  of  commands  and  controls  is  equal,  and  the  compound 

★ 

matrix  in  Eq.  16  is  nonsingular,  then  it  can  be  inverted,  and  Lx  and 
★ 

Au  are  defined  by 
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Ax*  -  -  F’V-H  F^G  ♦  Hu)_1  Ay^ 
Au*  =  (-  H  F_1G  +  H  )_1  Ay, 
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(17) 

(18) 

(19) 


Throughout  this  report,  the  relations  in  Eqs.  18  and  19  will  be  written 
as 


* 


Ax 


S12 


‘S  '  S22 

k 


(20) 

(21) 


for  reasons  explained  in  Ref.  15. 


Substituting  Eqs.  20  and  21  into  Eq.  14  and  rearranging,  the 
control  inputs  (as  functions  of  the  pilot's  inputs)  are 


T 


^  “  S22  “  c(A4  -  Su\] 

A  block  diagram  of  this  system  is  shown  in  Fig .  2-1. 


(22) 


OUTPUT 

SELECTION 


Figure  2-1.  Type  0  CAS  Structure 


It  is  interesting  to  note  that  due  to  the  feedforward  of  y.  through 

k 

S>22»  i*  is  possible  to  have  zero  "hang-off"  error  even  though  this  is 
a  Type  0  system.  A  condition  of  zero  hang-off  depends,  however,  on  no 
disturbances  and  knowing  the  exact  system  parameters  when  calculating 
S22  and  S^2* 


The  gains,  C,  are  calculated  using  linear  quadratic  (LQ)  regulator 
theory  (Ref.  12).  The  objective  of  this  method  is  to  find  the  feedback 
gains  which  allow  the  control  law  to  minimize  the  quadratic  cost  functional. 


J 


(AxJ  QAac  +  AuT  RAu)  dt 


(23) 
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where  Q  is  positive  semi-definite  and  R  is  positive  definite.  Matrices 
Q  and  R  weight  the  importance  of  suppressing  perturbations  in  states 
and  control  inputs  away  from  the  commanded  values.  Therefore  Q  and  R 
are  design  parameters  which  provide  flexible  and  efficient  means  of 
changing  closed- loop  response  to  satisfy  the  response  criteria. 

Gains  calculated  to  minimize  Eq.  23  will  not  minimize  its 
sampled-data  counter-part , 

N-l  _  T 

J  =  Z  (Ax!  QAxt  ♦  Au!  RAu.)  (24) 

i=o  A 

for  the  same  Q  and  R.  Reference  12  outlines  a  method  of  picking  con¬ 
tinuous-time  weighting  matrices  that  minimize  the  continuous-time  cost 
functional  even  though  a  sampled-data  regulator  is  used.  In  this  case 
the  cost  functional  for  the  sampled-data  case  penalizes  the  system 
continuously  in  time,  i.e.,  not  just  at  the  sampling  instances.  Mini¬ 
mizing  a  continuous  cost  functional  with  a  sampled-data  regulator  is 
equivalent  to  making  the  state  trajectory  obtained  from  the  digital 
CAS  as  close  as  possible  to  that  obtained  with  a  continuous  controller. 
This  is  true  because  both  controllers  are  derived  from  the  same  per¬ 
formance  criterion,  J.  Using  the  method  of  Ref.  12,  Eq.  23  is  set 
equal  to  the  sampled-data  cost  functional: 

J  =  I  [Ax£  QAXj.  +  2AxJmAu^  ♦  Au^RAu^]  (25) 

k=c 

The  discrete  weighting  matrices  are  found  to  be: 


k+1  *TCT)  Q  $(T)  dt 


/  k+1  T 

M  =/  1  4>  (T)  Q  T(T)  dt 


r  =J  k+1[R  +  rT 


(T)  Q  T  (T)  dt 


The  matrices  Q  and  R  are,  respectively,  positive  semi-definite  and 
positive  definite.  The  steady-state  discrete  Riccati  equation  is, 


p  =  $tp$  _  (rTp$+MjT(R+rTpr)"1(rTp$+M)  +  q 


where  P  is  the  equation's  steady-state  solution.  The  gains  C,  are 


given  by: 


c  =  (R  +  rTpr)"1(rTp$  +  m) 


To  find  the  gain  matrix  C  for  a  certain  Q  and  R,  Q,  M  and  R  must  be 
calculated  using  Eqs,  26-28.  Equation  29  then  is  solved,  and 
these  results  are  used  in  the  gain  equation  (Eq.  30)  to  compute  C . 


In  summary,  this  section  has  presented  a  way  to  control  a 
continuous-time  system  with  a  discrete-time  Type  0  CAS,  The  weighting 
matrices  Q  and  R  can  be  specified  for  a  continuous-time  cost  functional 
but  implemented  with  a  discrete-time  control  law  (Eq.  22).  This  is 

A  A  A 

easier  than  trying  to  specify  Q,  M,  and  R  directly,  because  Q  and  R 


can  be  based  on  state-control  tradeoffs  in  the  continuous-time 
problem.  Also,  Q  and  R  can  be  picked  as  diagonal  matrices  whereas 

A  A 

Q  and  R  would  not  normally  be  diagonal. 

2.2  EFFECTS  OF  THE  SAMPLING  INTERVAL 

As  mentioned  in  Chapter  1,  a  digitized  analog  control  law 
normally  requires  very  high  sampling  rates  which  could  put  excessive 
computational  burden  on  a  microprocessor-based  system.  High  sampling 
rates  are  necessary  for  mimicked  analog  systems  because  of  the  phase 
lags  associated  with  zero-order  hold  sampling.  Reference  16  pre¬ 
sents  a  graph  (Fig.  2-2)  that  shows  the  relation  between  phase  lag 
and  sampling  frequency  when  the  signal  is  first  passed  through  a 
second-order  anti-aliasing  filter  with  20  dB  attenuation  at  the 
sampling  frequency,  fsample. 


Figure  2-2.  Frequency  Ratio  Versus  Phase  Lag  for  20  dB  Attenuator. 
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For  instance,  if  a  phase  shift  of  10  deg  is  tolerable,  the  sampling 
frequency  should  be  approximately  30  times  the  highest  frequency  of 
interest,  according  to  Fig.  2-2. 


The  constraints  imposed  by  the  phase  lag  problem  are  eliminated 
when  the  sampled-data  regulator  theory  of  Section  2.1  is  used  in  the 
controller's  design.  It  is  not  too  difficult  to  see  that  sampling 
lag  effects  are  taken  into  account  with  a  sampled-data  regulator  because 
the  cost  functional  is  being  minimized  in  continuous-time  and  not  just 
at  discrete  instances.  Not  quite  so  obvious  is  why  it  is  not  necessary 
to  sample  at  higher  than  twice  the  highest  frequency  mode  to  be  con¬ 
trolled.  This  constraint  is  avoided  by  taking  advantage  of  our  prior 
knowledge  of  the  system  in  developing  the  control  law.  As  can  be  seen 
from  Eqs.  22-30,  the  aircraft's  dynamics  and  control  effects  are  in¬ 
cluded  right  in  the  control  law.  For  any  given  control  input,  the 
digital  CAS  "knows"  what  the  aircraft's  response  will  be  and  controls 
it  to  match  y,  to  the  degree  specified  by  the  weighting  matrices. 


There  is,  of  course,  a  maximum  value  for  the  sampling  inter¬ 
val  for  any  given  system  to  be  controlled.  This  maximum  can  be  dictated 
by  one  of  two  different  constraints.  The  two  constraints  are  1)  error 
build-up  between  samples,  and  2)  pilot  acceptance.  These  two  limiting 
factors  are  explained  below. 
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Error  Build  Up  Between  Samples 

If  there  were  no  disturbance  inputs  to  the  system  and  the 
dynamic  and  control  parameters  of  the  system  were  known  exactly,  then 
this  constraint  would  not  be  present.  But  in  the  real  world  there  are 
disturbances  and  parameter  variations,  so  errors  between  the  desired 
and  actual  states  of  the  aircraft  can  build  up  while  the  system  is 
running  "open-loop"  between  the  discrete  control  inputs.  The  sampling 
interval  should,  therefore,  be  selected  so  that  error  build  up  is 
limited  to  acceptable  values  under  "worst  case"  conditions  (Ref.  17). 

The  error  buildup  can  be  described  by  propagating  the  state  covariance 
with  the  differential  equation, 

P  =  FP  +  PFT  ♦  W  (31) 

where  P  is  the  state  covariance  matrix  and  W  is  the  disturbance  covari¬ 
ance  matrix  (Ref.  18). 

Error  build-up  because  of  too  long  a  sampling  interval  can  lead 
to  large  control  change  increments,  which  may  be  damaging  to  the  control 
surface  actuators  and  associated  mechanisms.  These  large  changes  are 
commanded  instantaneously,  which  causes  the  actuators  to  deflect  the 
controls  at  the  maximum  rate.  These  quick,  pulsing  movements  may  set 
up  vibrations  and  high  stresses  on  the  actuators  and  linkage.  This 
may  effectively  reduce  the  operating  life  of  these  components. 

Pilot  Acceptance 


Because  a  DFCS  changes  its  output  only  at  discrete  instants. 


the  plane’s  motions  set  by  a  pilot's  commands  may  not  occur  until  up 
to  one  sampling  interval  after  the  command.  If  the  pilot  moves  a 
control  and  the  aircraft's  response  is  delayed  perceptably,  his 
ability  to  control  is  degraded.  The  sampling  rate  must  be  at  least 
high  enough  so  the  pilot  perceives  no  appreciable  time  delay  in 
response  to  his  inputs,  as  defined  for  each  phase  of  flight.  It  is 
necessary  to  conducted  piloted  experiments,  either  in  flight  or  in 
ground-based  simulators,  to  determine  this  significance  of  time-delay 
effects . 

2.3  OBJECTIVES  OF  CONTROL  DESIGN 

Two  different  design  objectives  for  the  control  law  are 
addressed,  both  aimed  at  modifying  the  basic  short-period  response 
of  the  aircraft  in  straight-and-level  flight.  The  first  is  to  improve 
the  response  of  the  commanded  state  to  a  step  input  from  the  pilot. 

An  improved  response  is  judged  by  a  decrease  in  rise  time  and  over¬ 
shoot  (an  increase  in  the  longitudinal,  closed-loop,  short-period 
natural  frequency  and  damping).  The  degree  to  which  the  response  can 
be  improved  is  constrained  by  the  maximum  rate  of  deflection  of  the 
control  surface  and  the  requirement  to  satisfy  "Level  1"  flying  quality 
specifications  for  the  short-period,  as  defined  by  MIL-F-8785B(ASC.) 
(Ref.  19).  The  FBW  system  in  the  VRA  is  set  to  disengage  automatically 
if  the  error  between  actual  and  commanded  deflections  of  elevator  is 
greater  than  approximately  5  deg.  The  error  to  be  tolerated  in  the 
analytic  study  is  taken  to  be  this  value. 
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The  specifications  for  the  short-period  motion  in  the  Military 
Specification  on  Flying  Qualities  are  stated  as  constraints  on  the 
natural  frequency  and  damping  of  that  mode.  For  cruise,  the  minimum 
damping  ratio  is  0.3  and  the  maximum  is  2.0.  Natural  frequency  limita¬ 
tions  are  expressed  in  a  graph  with  the  frequency,  ( u)n  )  related  to 

sp 

normal  acceleration  (nz)  per  angle  of  attack  (a)  (Fig.  2-3).  The  basic 
VRA  is  marked  on  Fig.  2-3  and  is  in  the  Level  1  flying  qualities 
area. 

It  is  not  always  desirable  to  speed  up  the  response  and  reduce 
overshoot  of  the  commanded  inputs.  Sometimes  it  is  necessary  to  slow 
the  response  down  because  the  basic  aircraft  is  too  quick  and  over¬ 
sensitive  to  commanded  inputs.  The  weighting  matrices  can  be  found 
to  do  this  by  a  method  of  model  following  explained  in  Section  2.6. 
Therefore,  the  second  design  objective  for  the  control  law  is  to  allow 
simulation  of  another  aircraft's  dynamics  by  matching  the  closed-loop 
short-period  modal  characteristics  of  the  VRA  with  the  aircraft  to  be 
modeled.  Success  of  the  match  is  checked  analytically  by  comparing 
stability  derivatives  of  the  closed-loop  and  model  F  matrices. 


In  this  case,  the  design  objective  for  augmenting  the  short-period  is 
not  to  increase  the  natural  frequency  or  damping  ratio,  but  to  show 
that,  in  fact,  the  short-period  poles  can  be  arbitrarily  placed  in  the 
left-half-plane  as  required  using  the  LQ  control  methods. 

2.4  PITCH  RATE  COMMAND  CONTROL  STRUCTURE 

Pitch  rate  (q)  is  a  desirable  state  for  longitudinal  stick 
(6s)  to  command  because  at  low  speeds,  when  performing  some  sort  of 
tracking  task,  q  is  the  motion  cue  the  pilot  senses  for  a  deflection 
of  <5s>  The  simple  Type  0  control  law  derived  in  Section  2.1  is  now 
applied  so  the  response  of  Aq  to  A6g  can  be  modified  as  required. 

The  procedure  begins  by  defining  the  dynamics  that  give  changes  in  q. 

No  coupling  between  the  longitudinal  and  lateral-directional 
dynamics  will  be  considered,  so  the  state  and  control  vectors  of 
Eqs.  2  and  3  reduce  to: 

xT  =  [  u  w  q  8  ]  (32) 

uT  =  [  6e  6t  6p  ]  (33) 

States  u  and  w  are  not  directly  measurable  by  the  sensors  of  the  VRA, 
so  x  is  redefined  as, 

xT  =  [  V  Y  q  a  ]  (34) 


where  the  velocity,  V,  and  climb  angle,  y  are  taken  from  the  velocity- 
axis,  and  a  is  taken  from  the  wind-axes.  The  F  matrix  associated  with 


the  above  state  vector  is  derived  from  the  nonlinear  equations  of 
motion  (Ref.  13)  as  in  Eq.  8  and  is, 


F  = 


v 

“v 


-g  cos  Yo  Dq 


g/Vosin  Yc 


M. 


-g/Vosin  yo  (1-Lq) 


-L 


a 


(35) 


where  the  elements  of  F  are  stability  derivatives  and  inertial  effects 
that  represent  the  sensitivity  of  state  rates  to  state  perturbations. 
At  105  KIAS  in  straight-and- level  flight,  the  VRA's  F  matrix  is: 


F  = 


-.073 

.002 

.002 

-.002 


-32.2 

0 

0 

0 


0 

0 

■2.08 

1 


0 

2 

-8.35 

-2 


(36) 


where  angles  are  in  radians,  velocity  in  fps,  and  acceleration  is  in 
2 

fps  .  The  F  matrix  can  be  divided  into  four  (2  x  2)  matrices  that 
contain  the  parameters  for  the  short-period  and  phugoid  modes: 


F  = 


Phugoid 

Parameters 


Phugoid  to 

Short-Period 

Coupling 


Short-Period 

to 

Phugoid  Coupling 


Short-Period 

Parameters 


(37) 


Since  the  off-diagonal  matrices  have  elements  with  values  of  near  zero, 
the  short-period  is  described  mostly  by  changes  in  q  and  a,  and  the 
phugoid  is  described  mostly  by  changes  in  V  and  Y.  This  implies  control 
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of  q  response  can  be  accomplished  by  examining  a  reduced  second-order 
model  of  the  VRA,  using  the  (2  x  2)  submatrix  of  Eq.  37  with  the  short  - 
period  parameters.  The  controls  vector  is  reduced  to  <5g  in  order  to 
simplify  the  initial  control  law  design.  The  G  matrix  is  (2  x  1)  and 
has  the  elements, 


for  the  same  flight  condition  as  Eq.  36.  The  resulting  model  used  for 
the  design  of  the  pitch  rate  controller  is: 


AZ  =  Hx  Ax  ♦  Hu  A6e 


(40) 


The  reduced-order  VRA  dynamics  along  with  a  choice  of  weighting 
matrices,  Q  and  R,  can  be  used  to  generate  the  gains  C  and  matrices  Sj2 
and  S22.  Appendix  B  contains  an  explanation  of  the  APL  functions  used 
to  find  and  evaluate  C,  Sj2  and  S22.  Essentially  the  procedure  is  as 
shown  in  Fig.  2-4.  First,  matrices  H*  and  H„  are  defined  to  specify 
what  variables  the  pilot  will  command  with  in  this  case: 


Hx  =  [1  0]  (41) 

Hu  =  0  (42) 
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Figure  2-4.  Sequence  for  Finding  Optimal  Gains. 
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The  sampling  rate  is  chosen  next  because  it  is  used  in  Eqs.  11  and  12 
to  find  4>  and  T.  and  S22  are  calculated  from  Eqs.  18  and  19  next 
since  they  do  not  depend  on  how  the  states  and  controls  are  weighted. 
For  the  pitch  rate  CAS  these  matrices  are: 


S12 


’22 


.5 


Choosing  Q  and  R  will  set  the  closed- loop  response  of  the  system,  and 
changing  the  elements  of  these  matrices  allows  the  desired  response  to 
be  attained.  An  appropriate  starting  choice  for  the  weighting  factors 
in  Q  and  R  is  the  inverse  of  the  maximum  allowable  mean-square  values 
of  the  states  and  control  (Ref. 20): 


R-Va62  (43)  ,(44) 

E 

max 

In  the  iterative  technique  of  developing  the  desired  time  response, 
the  weighting  factors  can  be  changed  higher  or  lower  than  the  initial 

A  A  ^ 

set  as  necessary.  Matrices  Q,  R,and  M  are  found  by  performing  an  Euler 
integration  of  Eqs.  26-28.  The  Riccati  equation  (Eq.  29)  is  solved, 
yielding  the  matrix  P  necessary  to  derive  the  gains  from  Eq.  30.  The 
closed-loop  F  matrix,  F^.  is  the  continuous-time  equivalent  of 
and  can  be  calculated  from  the  expansion  of  (Ref. 21), 

fcl  •  T  >"  *CL  <45> 


‘V  0 
0  V 


which  is 


fCL=  f  -  1/2(»CL-I)2  *  1/3  C4>cl-I)3-  ...] 


(46) 


By  substituting  Eq.  14  into  Eq.  10,  4>£l  can  easily  be  found  to  be: 


$CL  =  *  -  r  c 


The  eigenvalues  of  the  closed- loop  system  are  the  roots  of  the  charac¬ 
teristic  equation  specified  by  the  determinant: 


|sl  -  fcl|  =  0 


Continuing  the  procedure,  time  histories  are  generated  using  Eqs.  10 
and  14.  If  examination  of  the  time  response  and  closed-loop  eigen¬ 
values  does  not  show  the  desired  response,  it  is  necessary  to  go  back 
and  re-choose  the  matrices  Q  and  R. 

As  mentioned  previously  in  Section  2.3,  there  is  no  established 
general  relation  between  quadratic  weightings  and  time  response  for  any 
given  controller  or  plant.  The  particular  relation  between  response  to 
a  step  input  (measured  in  rise  time  to  95  percent  of  commanded  value, 
and  percent  overshoot  of  commanded  value)  and  the  Q  matrix  weightings 
is  shown  in  Fig.  2-5.  Matrix  R  is  held  constant  with  a  value  of 
l/(.17rad)  =  32.8,  and  the  sampling  rate  is  10  sps.  The  curves  radi¬ 
ating  outward  correspond  to  decreasing  rise  times  and  increasing  initial 
control  deflections.  The  straight  lines  intersecting  the  origin  repre¬ 
sent  decreasing  overshoots  when  approaching  the  Aq  axis.  The  dotted 
line  represents  the  approximate  boundary  between  a  closed-loop  system 
with  complex  eigenvalues  (below  the  line)  and  one  with  real  eigenvalues. 
About  25  data  points  were  used  to  derive  the  curves  of  Fig.  2-5  using 
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Figure  2-5.  Relation  Between  Weighting  and 
Time  Response  -  q  CAS;  R  =  32.8. 


interpolation.  The  general  relationship  is  shown,  and  if  a  particular 
response  is  desired  of  this  CAS,  Fig.  2-5  can  give  very  good  starting 
values  for  the  weightings,  requiring  only  a  few  iterations  to  find 
the  necessary  Q.  The  relation  between  weightings  and  response  does 
not  vary  much  at  all  for  different  sampling  rates  because  sampled-data 
control  methods  minimize  the  cost  functional  continuously  in  time 
regardless  of  the  sampling  interval.  Figure  2-5  can  therefore  be  used 
for  sampling  rates  other  than  10  sps. 

Time  history  plots  for  a  few  representative  weightings  were 
generated  on  Princeton  University's  IBM  370/158  time=sharing  computer 


nit 


) 


(Figs.  2-6  through  2-11).  All  plots  are  for  a  step  command  input  of 
5  deg/sec  using  a  full  fourth-order  model  of  the  VRA's  longitudinal 
dynamics.  The  Aq  response  on  the  open- loop  aircraft  (Fig.  2-6)  has 
a  47.5%  overshoot  and  0.19  sec  rise  time.  As  the  velocity  decreases, 

Aq  falls  off  because  there  is  no  throttle  input.  Angle  of  attack 
shows  the  expected  initial  increase  as  the  phugoid  oscillation  sets  in. 
Weighting  factors  in  the  Q  matrix  of  (25,25)  show  a  definite  improve¬ 
ment  in  rise  time  and  overshoot  of  the  Aq  response  (Fig.  2-7).  Because 
the  CAS  is  Type  0,  and  there  is  no  feedback  to  throttle,  Aq  again  shows 
a  decrease  but  not  quite  as  much  as  the  open- loop  case.  The  elevator 
initially  attains  a  larger  than  steady-state  deflection  to  improve 
rise  time,  but  then  dips  below  steady-state  value  to  control  the  over¬ 
shoot.  Finally,  the  elevator  assumes  the  steady-state  deflection  for 
Aq  *  5  deg/sec.  The  plots  for  weightings  of  (50,50)  and  (75,25) 

(Figs.  2-8  and  2-9)  show  further  improvements  in  time  response,  but 
the  initial  elevator  deflection  is  larger.  The  time  responses  with 
weightings  (25,25)  and  (75,25)  are  also  plotted  (Figs.  2-10  and  2-11) 
for  a  sampling  rate  of  20  sps.  It  can  be  seen  that  they  have  the 
same  state  trajectories  as  the  10  sps  plots  of  the  same  weightings. 

The  control  trajectories  are  different  because  the  CAS  is  taking  into 
account  the  difference  in  phase  lag  between  the  two  sampling  rates. 

The  gains,  resulting  step  responses,  and  closed-loop  eigen¬ 
values  obtained  from  the  simulation  are  given  in  Table  2-1.  A- nega¬ 
tive  gain  implies  negative  feedback  while  positive  gains  give  positive 
feedback.  The  phugoid  eigenvalues  change  very  little  for  different 
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Figure  2-8.  Pitch  Rate  Command  Step  Response  for  Q  Weightings  50,50 
(yd  =  5°/ sec,  10  sps,  105  K1AS) . 
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re  2-10.  Pitch  Rate  Command  Step  Response  for  Q  Weightings  25,25 
(yd  =  5°/sec,  20  sps,  105  KIAS) . 
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Figure  2-11. 
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Pitch  Rate  Command  Step  Response  for  Q  Weightings  75,25 
(yd  «  5°/sec,  20  sps,  105  KIAS) . 
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TABLE  2-1 


Gains,  Predicted  Step-Response  and  Eigenvalues  for  q  CAS  (105  KIAS) 


Open  Loop 


10  sps 


20  sps 


Q 

Weightings 

Gains 
q  a 

Eigenvalues 

Short-Period, 

Phugoid 

Response 

Rise  Time,  Overshoot, 
Sec  % 

- 

- 

-2. 04+j  2. 89 
-.033+j .228 

0.19 

47.5 

[25,25] 

-.427 

.324 

-10.74,-2.68 

-0.45*j.l39 

0.15 

14.7 

[50,50] 

-.554 

.360 

-16.96,-2.43 
-.0^-6+j  .123 

0.09 

11.8 

[75,25] 

-.620 

.468 

-23.24,-2.20 

-.048+j.l04 

0.08 

6.7 

[25,25] 

-.541 

.283 

-10.30,-2.68 

-.045+j.l33 

0.14 

14.3 

[75,25] 

-.876 

.439 

-19.31,-2.20 
- .049+p .092 

0.1 

6.4 

gains  due  to  the  weak  coupling  between  short-period  and  phugoid  modes. 
The  un  is  increased  from  the  open- loop  value  of  3.54rad/sec  (Table 
2-2)  for  the  3  weightings  examined.  This  increase  is  primarily  seen 
in  the  decrease  in  rise  time  for  the  commanded  state,  Aq,  and  is  con¬ 
sidered  an  improvement  by  the  time  response  criterion.  However,  re¬ 
ferring  to  Fig.  2-3,  the  closed- loop  for  weightings  of  (50,50) 

sp 

and  (75,75)  is  outside  the  Level  1  area  for  n/a  constant  at  9.3g/rad. 
This  implies  the  flying  qualities  of  the  VRA  may  be  unsatisfactory  in 
actual  flight  tests  if  gains  calculated  from  these  weightings  are  used. 


» 
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TABLE  2-2 

Closed-Loop  and  CSp  for  q  CAS  (105  KIAS) 


Q 

Weightings 

^sp 

(rad/ sec) 

^sp 

Open  Loop 

3.54 

0.576 

[25,25] 

5.39 

1.25 

[50,50] 

6.45 

1.51 

[75,25] 

7.19 

1.77 

Examination  of  F^l,  calculated  from  Eq.  46,  will  show  how 
the  stability  derivatives  have  been  augmented  for  the  particular 
weightings  chosen.  For  example,  the  second-order  Fql  for  weightings 

of  (75,  25)  is; 


M„ 

M„ 

-23.5 

-4.48 

a 

S 

1-Lq 

-L 

a 

1.40 

-1.93 

- 

- 

Comparing  this  to  the  open- loop  F  matrix  of  Eq.  39,  M  has  increased 
by  approximately  a  factor  of  10,  has  decreased  by  a  factor  of  two, 
and  some  negative  lift  due  to  q  has  developed.  The  stability  deriva¬ 


tive,  Lq  has  not  changed  much  because  the  only  control  surface  used 
by  the  CAS  is  elevator  and  L^  is  very  small.  To  augment  Lq,  direct 
lift  control  from  the  flaps  is  needed.  The  Fql  for  weightings  of 


(25,25)  and  (50,50)  are, 


-11.5  -5.93 


1.09  -1.96 


•17.5  -6.56 


1.22  -1.92 


respectively. 


2.5  NORMAL  ACCELERATION  COMMAND  CONTROL  STRUCTURE 


Normal  acceleration  commanded  by  6s  is  the  next  CAS  structure 
examined.  The  development  is  very  similar  to  that  of  the  pitch  rate 
CAS  of  Section  2.4;  hence,  the  principal  differences  are  given  here. 
The  gains  are  again  calculated  with  a  second-order  model  but  using 
the  states  q  and  nz.  Using  the  relation, 


n  =  Z  a  ♦  Zt  6 
z  a  Oj;  E 


the  model  represented  by  Eq.  39  can  be  used  if  the  correct  substitu¬ 
tion  for  Act  is  made.  Since  the  effects  of  Z$g  are  so  small,  the 
second  term  of  Eq.  50  is  dropped  to  make  the  substitution  easier. 
Making  the  change,  the  second-order  F  matrix  is: 


VLa  -La 


-2.08  -.0235 


354.9 


Normalizing  acceleration  by  the  force  of  gravity,  the  resulting  model 


» 
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is  given  as; 
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Figure  2-12.  Relation  Between  Weightings  and 
Time  Response  -  nz  CAS;  R  =  32.8. 

Time  history  plots  for  a  few  representative  weightings  are 
presented  in  Figs.  2-13  through  2-16.  All  plots  are  generated  with 
a  step  command  input  of  0.25g  using  a  full  fourth-order  model  of  the 
VRA's  longitudinal  dynamics.  The  Anz  response  for  the  open- loop  air¬ 
craft  (Fig-  2-13)  has  a  10.8%  overshoot  and  0.7  second  rise  time. 
Velocity  is  falling  off  and  holds  constant  at  1.3  deg  because  its 
position  is  only  a  function  of  the  step  input  on  Weightings  of 

(5,10)  in  the  Q  matrix  and  R  =  33  (Fig.  2-14)  shows  a  large  improve¬ 
ment  in  rise  time  for  Anz  response.  However,  now  Aq  overshoots  its 
steady- state  value  quite  a  bit  more  than  the  open- loop  case.  The 
elevator  initially  deflects  more  than  5  deg,  which  indicates  the  FBW 


TIME  (sec) 


TIME  (sec) 


system  of  the  VRA  will  disengage  for  a  .25g  true  step  input. 

Figure  2-15  shows  the  response  for  a  heavier  weighting  (R=250)  on 
Au.  The  initial  step  of  A6g  is  now  less  than  3  deg  but  rise  time 
on  A"z  has  been  slowed  by  0.2  sec.  The  peak  on  the  Aq  overshoot 
also  is  less  with  the  heavier  control  weighting.  A  set  of  plots 
drawn  for  the  same  weightings  as  Fig.  2-14,  but  a  sampling  rate  of 
20  sps  (Fig.  2-16),  shows  once  again  that  a  change  in  sampling  rate 
does  not  effect  the  state  trajectories  if  new  gains  are  calculated 
to  take  into  account  the  change  in  phase  lag  of  the  inputs  to  the 
digital  CAS.  The  control  trajectory  is  different  from  Fig.  2-14  as 
a  result. 

The  gains,  resulting  step  responses,  and  closed-loop  eigen¬ 
values  obtained  from  the  computer  simulation  are  given  in  Table  2-3. 

TABLE  2-3 


Gains,  Predicted  Step-Response  and  Eigenvalues  for  nz  CAS  (105  KIAS) 


Weightings 

Q. 

R 

Gains 

q 

n_z 

Eigenvalues 
Short  Period, 
Phugoid 

Response 

Rise  Time,  Overshoot, 
Sec  % 

Open  Loop 

-  2. 04+ j 2. 89 
-.033+;. 228 

0.7 

10.8 

[5,10] 

33 

-.533  -. 

181 

-6.47+j6.04 
- . 036+j . 252 

0.35 

3.0 

[5,10] 

250 

-.238  -. 

0535 

-3.71+j4.09 
-. 0360+ j. 257 

0.52 

5.6 

[5,10] 

33 

-.620  -. 

250 

-6 . 52+ j  5.95 
-,0364+j .252 

0.35 

3.15 
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Figure  2-16.  Normal  Acceleration  Step  Response  for  Q  Weightings  5,10 
R  *  33  (yd*.25g,  20  sps ,  105  KIAS) . 
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The  au  for  the  weightings  with  R  =  32.8  and  250  is  8.85  rad/sec 
sp 

and  5.52  rad/sec  respectively.  This  increase  in  ianSp  from  open- 
loop  (3.54  rad/sec)  accounts  for  the  decrease  in  rise  time  for 
Anz.  This  is  considered  an  improvement  by  the  time  response  criterion 
but  for  the  case  where  t%Sp  =  8.35,  it  no  longer  meets  the  require¬ 
ment  set  by  the  Military  Specification  on  Flying  Qualitites  (Ref. 20). 
Referring  to  Fig.  2-3,  the  augmented  Un  puts  the  aircraft  in  the 

bp 

Level  2  area  of  the  chart  for  n/c*  still  equal  to  9.3  g/rad.  This 
indicates  the  control  law,  with  gains  calculated  from  weightings  of 
and  R  =  32.8,  may  not  be  satisfactory  to  the  pilot  in  the 
actual  flight  tests. 

The  closed-loop  stability  derivatives  may  be  examined  by 
calculating  FCL  from  Eq.  46.  The  Fq^  for  the  systems  with  R  *  32.8 
and  R  =  250  are 


-11.5 

-5.33 

and 

’-5.49 

- 1 . 8o" 

11.53 

-1.48 

11.09 

-1.93 

_ 

- 

respectively.  These  can  be  compared  to  the  open-loop  F  matrix  of 
Eq.  52  to  see  how  the  stability  derivatives  have  changed. 

2.6  DYNAMIC  MODEL  MATCHING 

The  method  of  choosing  weighting  matrices  Q  and  R  described 
in  the  preceding  two  sections  can  result  only  in  decreases  in  over¬ 
shoot  and  rise  time  for  the  time  response  of  the  commanded  state; 
however,  it  may  be  desirable  to  change  the  modal  characteristics  of 
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an  aircraft  in  the  opposite  way  if  the  basic  plane's  response  is  too 
fast  or  too  well  damped.  This  can  be  done  with  the  previously  des¬ 
cribed  LQ  control  methods  by  using,  in  addition  ,  the  method  of 
implicit  model- following. 

The  procedure  begins  by  defining  the  eigenvalues  desired 
for  the  closed- loop  system.  Next,  the  F  matrix  associated  with  the 
eigenvalues  through  Eq.  48  is  found.  By  the  method  of  implicit  model¬ 
following,  it  is  possible  to  find  the  weightings  which  give  the  closed- 
loop  system  the  correct  F  matrix  (Ref.  22).  Just  as  in  Section  2-1 
the  design  is  posed  as  a  linear-quadratic  control  problem,  but  with  a 
different  quadratic  cost  functional 

CD 

J  =  f  (AxTQAx  +  2  A£TMAu  +  AuTRAu)  dt  (56) 

■•o 

where 

Q  *  (F-FM)TW  (F-Fm)  (57) 

M  «  (F-Fm)T  W  G  (58) 

R  «  GT  W  G  (59) 

The  matrix  W  is  a  diagonal  weighting  matrix  and  F^  is  the  F  matrix 
to  be  modeled.  If  the  modeling  was  to  be  done  in  continuous-time,  the 
Riccati  equation  would  be  solved  next  and  the  feedback  gains  calcu¬ 
lated.  However,  in  the  discrete-time  case,  the  weightings  must  be  modi¬ 
fied  to  account  for  the  effective  lag  of  the  sampling  process.  Unlike 
the  cost  functional  of  Eq.  23,  there  is  a  cross  weighting,  M,  in  Eq.  56 


which  cannot  be  transformed  to  its  discrete  counterpart  just  by  using 
Eqs.  26-28  of  the  procedure  outlined  in  Section  2.1.  Following  the 
development  in  Ref.  IS,  a  transformation  from  the  continuous -time 
weighting  matrices  (Q,  M,  and  R)  to  the  discrete-time  matrices  (Q,  M, 
and  R)  was  derived  (see  Appendix  A).  The  discrete  weighting  matrices 
are  found  to  be 


*  ftk+l  T 
Q  =  I  **  Q  4>  dt 

'*k 

„  (tk+l 

M  =  (rT  Q<J>  +  M*)dt 


^  ftk+l  T 

R  =  (R  +  Tt  QT  +  2  MT)dt 


The  gains  are  calculated  (Eq.  30)  from  the  solution  of  the  Riccati 
equation  (Eq.  29),  and  is  found  with  Eq.  44.  If  Fql  is  not  exactly 
as  desired,  then  matrix  W  should  be  changed  and  a  new  set  of  quadratic 
weightings  calculated  (Fig.  2-17).  If  changing  W  does  not  yield  sat¬ 
isfactory  results,  the  F^  used  in  the  calculation  can  be  changed. 

APL  functions  to  perform  these  calculations  are  listed  and  explained 
in  Appendix  B. 

For  example,  say  it  is  necessary  to  find  the  feedback  gains 
for  a  closed- loop  system  with  u>nSp  «  1.56  and  CSp  =  0.76.  The  corre¬ 
sponding  eigenvalues  are  -1.18  ^  j  1.02.  The  FM  matrix  which  gives 
these  eigenvalues  must  be  found  with  the  constraint  that  some  elements 
of  F  for  the  VRA  cannot  be  changed  much  with  only  elevator  control. 
Specifically,  Lq  and  Lq  cannot  be  augmented  with  only  6E,  so  the  top 


elements  of  F  (Mq  and  M  )  are  the  values  to  be  found.  Using  Eq.  48, 
is  found  to  be: 


rM 


Choosing  the  diagonal  weighting  matrix  W  to  initially  be  a  (2x2)  identity 
matrix,  I_,  and  completing  the  rest  of  the  method  outlined  in  Fig.  2-17, 


fCL  = 


-.375 

.  1.00 


-1.68 

-2.01, 


For  this  particular  example, it  is  not  possible  to  make  F^g  close  to 
Fm  just  by  changing  W,  so  the  F^  used  in  the  calculations  is  changed 
instead,  keeping  W  *  I-.  After  a  few  iterations,  it  is  found  that, 


FM  = 


will  yield, 


CL  " 


For  reference,  the  weighting  matrices  and  resulting  gains  are  listed 
in  Table  2-4. 


The  next  chapter  outlines  the  hardware  and  software  aspects 
of  the  Micro-DFCS.  It  describes  how  the  control  laws  derived  in 
this  chapter  are  actually  implemented,  and  gives  test  results  and 
analysis  for  a  hybrid  simulation. 
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TABLE  2-4 

Weighting  Matrices  and  Gains  for  Implicit  Model-Following  Example 


Continuous-Time 

Weightings 

Discrete-Time 
Weightings  and  Gains 

r3.03  11. Si 

„  f-338 

.9651 

Q  = 

Q  = 

L 11-5  43.7J 

L.  965 

2. 83  J 

M  =  [21.8  82.6] 

M  =  [2.08 

6.03] 

R  =  156.3 

R  =  12.9 

C  =  [.161 

.467] 

3. 


CONTROL  LAW  IMPLEMENTATION  AND  HYBRID  TESTING 


The  control  law  developed  in  Chapter  2  is  coded  in  the 
software  of  the  Micro-DFCS  for  implementation  on-board  the  VRA. 

The  type  and  configuration  of  the  Micro-DFCS  equipment  specifies, 
to  a  large  extent,  in  what  manner  and  how  well  these  control  laws 
will  operate  with  the  existing  VRA  FBW  system.  The  efficiency  and 
versatility  of  the  equipment  and  software  determines  the  complexity 
of  control  laws  that  may  be  used.  The  control  law  designs  therefore 
have  an  impact  on  the  equipment  selected  and  organization  of  the  soft¬ 
ware.  This  chapter  outlines  some  of  the  specific  aspects  of  the  Micro- 
DFCS  that  support  the  implementation  of  the  control  design.  The  equip¬ 
ment  and  software  are  tested  with  an  analog  simulation  of  the  VRA  as 
a  logical  step  before  actual  flight  tests. 

3.1  EQUIPMENT  CONFIGURATION  OF  THE  MICRO-DFCS 

The  primary  task  of  the  Micro-DFCS  hardware  is  to  accept 
analog  information  to  calculate  the  control  laws,  and  to  generate 
analog  commands  for  aircraft  controls  at  periodic  instants  in  time. 

This  definition  of  the  computer's  main  function  suggests  some  uesirable 
characteristics  for  the  hardware  to  possess.  Since  the  control  law 
must  be  calculated  in  short  intervals  of  time,  the  computations  must 
be  very  fast  but  also  have  precision  at  least  as  great  as  the  reso¬ 
lution  of  the  digital  data  words  generated  by  the  analog  to  digital 
(A/D)  converter.  The  computational  speed  is  determined  by  the 
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instruction  cycle  time,  the  bit  length  of  the  data  word,  and  the 
speed  of  the  mathematical  operations.  The  A/I)  and  digital 
to  analog  (D/A)  converters  should  have  as  much  resolution  as 
the  sensors  will  provide  and  the  control  surfaces  will  move  (real 
sensors  and  actuators  do  not  have  infinite  precision  because  of  stic- 
tion).  An  interrupt  structure  with  at  least  one  interrupt  which  can  be 
.triggered  by  a  resetable  timer  is  needed  to  initiate  the  control  law 
calculation  at  precise  instances  in  time.  Other  factors  to  consider 
include  the  cost,  the  range  of  compatible  components  which  may  be 
used  to  expand  or  upgrade  the  computer  at  a  latter  date,  the  standard¬ 
ization  and  flexibility  of  the  bus  structure  (which  determines  the 
availability  of  components  compatible  with  the  system),  the  power  re¬ 
quirements,  and  the  software  available  to  support  the  microcomputer. 

The  Micro-DFCS  for  this  study  is  based  on  the  Intel  Single 
Board  Computer  (SBC)  using  the  Intel  8085  Central  Processing  Unit  (CPU) 
(Ref.  25).  The  SBC  system  is  a  modular  design  where  the  main  functions 
of  the  computer  are  divided  onto  separate  circuit  boards.  These 
boards  may  be  combined  to  configure  the  desired  system.  Multiple 
CPU’s  may  be  used  in  a  system  when  more  processor  time  is  needed. 

This  is  not  the  case  at  this  time,  but  may  be  necessary  as  more  com¬ 
plex  command  modes  are  implemented.  As  shown  in  Fig.  3-1,  the  Micro- 
DFCS  is  built  around  the  SBC  80/05  central  processing  board.  This 
is  supported  by  a  high-speed  mathematics  unit,  random-access  and 
programmable  read-only  memory  (RAM  and  PROM),  A/D  and  D/A  conversion 
boards,  and  a  hand-held  control-display  unit  (CDU) .  The  CPU  works 
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SENSORS  CONTROLS 


Figure  3-1.  Micro-DFCS  Microcomputer  Configuration. 

with  80  different  machine  instructions,  has  an  8-bit  data  word  length, 
and  an  instruction  cycle  time  from  2  to  6.1  Msec  (depending  on  the 
instruction).  A  detailed  description  of  the  microcomputer  components 
is  contained  in  Appendix  D. 

The  Micro-DFCS  is  functionally  related  to  the  VRA's  FBW 
system  as  shown  in  Fig.  3-2.  The  Micro-DFCS  can  operate  in  parallel 
with  the  existing  analog  system  or  either  can  operate  separately. 

This  is  an  important  aspect  of  the  Micro-DFCS  installation  because, 
in  this  configuration,  it  is  possible  to  model  another  aircraft 
with  the  VRA  using  the  analog  system  and  then  test  a  digital  control 
design  (for  the  modeled  aircraft)  using  the  Micro-DFCS.  Two  pilots 
fly  each  test  flight  for  reasons  of  safety  and  experimental  effi¬ 
ciency.  The  safety  pilot  has  conventional  aircraft  controls  with 
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Figure  3-2.  Overview  of  the  VRA/Micro-DFCS  System. 


mechanical  linkage  connecting  to  the  control  surfaces.  The  Evaluation- 
pilot  operates  the  experimental  controls,  which  in  this  case,  include 
a  center  control  stick  and  the  CDU. 

The  VRA's  motion  is  measured  by  various  air  data  and  inertial 
sensors  and  sent  to  a  ground  station  via  42  channels  of  telemetry.  A 
more  detailed  explanation  of  the  VRA's  systems  and  ground  station  is 
contained  in  Appendix  C. 

3.2  THE  OPERATIONAL  MICROCOMPUTER  SOFTWARE 

The  calculations  for  the  control  laws  are  performed  by  a  set 
of  software  programs  that  execute  every  sampling  instant.  These  Flight 
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Control  routines  are  specific  to  their  function  in  order  to  minimize 
the  execution  duty  cycle.  The  Flight  Control  routines  are  supported 
by  two  sets  of  programs  called  the  Executive  and  Utility  routines. 

They  are  designed  to  be  very  flexible  allowing  new  Flight  Control 
routines  to  be  easily  added.  These  three  categories  together  make 
up  the  operational  software  entitled  CAS-1.  The  routines  of  CAS-1, 
arranged  in  Chapter  format, are  listed  with  memory  requirements  in 
Table  3-1.  The  Executive  routines  handle  initialization,  CDU  inter¬ 
face,  and  one  of  three  error  detection  methods  used  in  CAS-1.  This 
error  detection  routine  checks  the  contents  of  memory  every  50  samples 
to  insure  the  coding  of  CAS-1  has  not  changed.  A  second  error  detec¬ 
tion  method  checks  for  mathematical  errors  on  every  operation.  The 
third  method  indicates  when  a  Flight  Control  routine  is  executing  pro¬ 
perly  by  flashing  a  light  at  a  steady  rate  on  the  pilot's  panel.  A 
more  detailed  description  of  the  error  detectors  and  a  brief  descrip¬ 
tion  of  all  routines  is  given  in  Appendix  E. 

The  flexibility  of  the  Executive  and  Utility  routines  is 
derived  mainly  from  the  modular  construction  of  the  program.  Major 
and  often  performed  tasks  are  divided  into  subroutines  that  are  very 
general  in  nature.  All  user-definable  parameters  and  variables  that 
change  in  the  operation  of  CAS-1  are  stored  in  RAM  outside  the  body 
of  CAS-1.  Therefore  any  part  or  all  of  CAS-1  may  be  stored  in. PROM. 
Developing  new  Flight  Control  programs  for  the  Micro-DFCS  can  be  done 
quickly  and  efficiently  by  using  these  existing  routines.  A  detailed 
explanation  on  defining  new  keyboard  commands  and  flight  control 
modes  is  contained  in  Appendix  E. 
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TABLE  3-1 


CAS-1  PROGRAM  TABLE  OF  CONTENTS  AND  MEMORY  REQUIREMENTS 


|  Memory  Words 

1.  EXECUTIVE  ROUTINES 

1.1 

INITIALIZATION 

75 

1.2 

CDU  INTERFACE  AND  COMMAND  RECOGNITION 

125 

1.3 

MEMORY  CHECK 

80 

2.  UTILITY  ROUTINES 

2.1 

ANALOG  TO  DIGITAL  CONVERSION 

18 

2.2 

ENTRY  ERROR 

53 

2.3 

BLINK 

21 

2.4 

CLEAR  LINE 

8 

2.5 

CONSOLE  OUTPUT 

10 

2.6 

COUNT-UP  DISPLAY 

72 

2.7 

DECIMAL  TO  HEX  CONVERSION 

24 

2.8 

ERASE  BLOCK  MEMORY 

7 

2.9 

MATH  ERROR  PROCESSOR 

111 

2.10 

HEX  TO  DECIMAL  CONVERSION 

60 

2.11 

NUMERIC  INPUT 

66 

2.12 

INTERRUPT  COUNT 

19 

2.13 

LIMIT  ANALOG  OUTPUT 

45 

2.14 

MATH  UNIT  DRIVER 

38 

2.15 

MODE  CHANGE 

142 

2.16 

MOVE  4  WORDS 

12 

2.17 

SERIAL  OUTPUT 

12 

2.18 

CALIBRATED  STEP  INPUT 

294 

3.  FLIGHT  CONTROL  ROUTINES 

3.1 

DIRECT  CONTROL  SET  UP  MODE  19 

77 

3.2 

PITCH  RATE  CONTROL  SET  UP  MODE  0 

134 

3.3 

PITCH  RATE  CONTROL  SET  UP  MODE  1  ) 

3.4 

PITCH  RATE  CONTROL  SET  UP  MODE  2  f 

95 

3.5 

PITCH  RATE  CONTROL  SET  UP  MODE  3 ) 

3.6 

NORMAL  ACCELERATION  CONTROL  SET  UP  MODE  10  135 

3.7 

DIRECT  INTERRUPT  SERVICE  ROUTINE 

94 

3.8 

PITCH  RATE  INTERRUPT  SERVICE  ROUTINE 

323 

3.9 

NORMAL  ACCELERATION  INTERRUPT  SERVICE 

ROUTINE  388 

TOTAL  =2540 

The  selected  control  mode  (Pitch  Rate,  Normal  Acceleration, 
or  Direct)  determines  which  control  law  will  be  calculated  at  each 
sampling  instant  by  the  Flight  Control  routine.  The  control  law 
calculation  has  the  highest  priority  over  all  other  tasks  because 
the  sampling  interval  is  to  be  kept  constant.  Figure  3-3  presents  a 
diagram  of  this  priority  structure  and  a  description  of  the  duty  cycle 
for  a  nominal  sampling  interval  of  0.1  sec. 
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Figure  3-3.  CAS-1  Program  Organization 


The  control  law  calculation  operates  in  the  "foreground"  and  all  other 
tasks  are  performed  in  the  remaining  available  time  (in  the  "back¬ 
ground")  .  The  duty  cycle  shows  again  that  the  routine  containing  the 
control  logic  is  executed  exactly  every  0.1  second  and  any  other  rou¬ 
tines  fill  in  on  the  available  "wait  state"  time. 


The  components  of  the  software  development  system  (Fig.  3-4) 
allow  programs  to  be  entered  and  tested  efficiently.  The  ground  chasis 
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Flight  Research  Laboratory 


Figure  3-4.  Equipment  Layout  for  Micro-DFCS 
Software  Development. 


holds  the  microcomputer  boards  and  provides  them  with  power.  The 
routines,  when  first  coded  in  assembly  language,  can  be  entered  into 
the  microcomputer  through  the  keyboard-CRT  unit.  The  text  of  the 
routine  can  be  edited  to  a  limited  extent  on  the  microcomputer  before 
it  is  sent  to  Princeton  University's  IBM  370/158  computer  via  an 
acoustic  coupler  and  the  telephone  lines.  The  IBM  computer  provides 
more  sophisticated  text  editing,  cross-assembly  of  Micro-DFCS  code, 
and  permanent  storage  of  all  routines.  The  assembled  routines  can 
be  loaded  to  the  microcomputer  from  the  IBM  computer  and  then  debugged 
using  the  keyboard-CRT  unit.  A  further  description  of  the  software 
development  system  is  contained  in  Appendix  D. 


For  actual  flight  tests,  the  programs  are  transferred  from 
the  development  system  to  the  Flight  Control  Computer  Unit  (FCCU) 
mounted  in  the  VRA.  The  programs  can  be  transferred  either  on  PROM 
or  on  a  battery  back-up  4K  RAM  board.  The  battery  back-up  RAM  is 
the  most  convenient  method,  because  program  modifications  made  in 
RAM  are  much  easier  and  quicker  than  programming  a  new  PROM.  CAS-1 
presently  occupies  approximately  2.SK  words  of  memory.  When  CAS-1 
is  expanded  past  the  4K  limit,  the  routines  that  will  require  no 
changes  (the  Executive  and  most  of  the  Utility  routines)  should  be 
relocated  in  PROM.  The  Flight  Control  routines  will  probably  be 
modified  quite  often  and  should  be  kept  in  battery  back-up  RAM  if 
possible. 

3.3  HYRBRID  SIMULATION  TESTS  AND  RESULTS 

As  a  logical  step  between  the  analytic  development  and  the 
flight  tests,  the  Micro-DFCS  hardware,  software,  and  implemented 
control  laws  were  tested  by  simulating  the  VRA's  longitudinal  dyna¬ 
mics  on  an  EAI  TR-48  analog  computer.  The  voltage  levels  of  the 
signals  interfacing  the  analog  computer  to  the  A/D  and  D/A  channels 
of  the  microcomputer  were  scaled  to  be  of  the  same  magnitude  and 
sign  as  the  outputs  from  the  actual  sensors  in  the  VRA.  This  allows 
a  "real-time"  simulation  of  the  Micro-DFCS  which  insures  that  all 
scalings,  gains,  and  other  parameters  within  the  software  have  been 
set  correctly.  The  analog  computer  diagram  for  the  longitudinal 
dynamics  of  the  VRA  and  a  table  of  poteniometer  settings  is  given  in 
Appendix  C  (Fig.  C-5  and  Table  C-2), 
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The  simulation  used  the  full  fourth-order  longitudinal 
equations  of  motion  that  give  the  effects  of  both  the  short-period 
and  phugoid  modes.  The  simulated  VRA's  q,  a,  nz,  and  6S  signals 
were  connected  to  the  analog  input  channels  of  the  Micro-DFCS.  The 
signal  generated  by  the  Micro-DFCS  was  connected  to  the  control 
input  of  the  modeled  VRA.  The  states  of  the  VRA  and  the  6g  output 
from  the  Micro-DFCS  were  recorded  on  a  strip  chart  recorder. 

The  three  flight  control  modes  —  Direct,  Pitch  Rate,  and 
Normal  Acceleration  —  were  tested  with  different  feedback  gains 
and  sampling  rates  (called  control  configurations).  Each  was  in¬ 
vestigated  by  generating  a  step  input  on  6g  and  examining  the  response 
of  the  state  variables.  Specifically,  the  rise  time  and  overshoot 
of  the  commanded  state  were  compared  to  those  predicted  by  the  ana¬ 
lytical  results.  The  step  of  6s  was  generated  internally  in  the  micro¬ 
computer  using  the  calibrated  step  routine  explained  in  Appendix  E. 

This  insures  that  a  true  step  is  made  for  each  input,  and  it  elimi¬ 
nates  differences  in  resulting  plots  that  could  be  due  to  the  input 
of  6S. 

The  control  law  calculations  were  timed  during  these  tests, 
and  it  was  found  that  the  Direct  Mode  routine  takes  S  msec  to  exe¬ 
cute,  the  Pitch  Rate  Mode  takes  19  msec,  and  the  Norma]  Acceleration 
Mode  takes  6  msec.  The  Normal  Acceleration  Mode  routine  executes 
faster  than  the  Pitch  Rate  Mode  because  the  former  is  coded  more 
efficiently  to  save  execution  time  and  not  memory  space.  This  is 
explained  in  more  detail  in  Appendix  E. 
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There  are  time  plots  at  the  end  of  this  section  for  15 


different  variations  of  the  three  basic  flight  modes.  Each  column 
of  Fig.  3-5  represents  one  run.  The  flight  control  configuration, 
and  the  Ay^  used  for  each  run  are  labeled  at  the  top  of  each  column. 

The  different  control  modes  are  listed  in  Table  3-2.  Each  mode  is 
for  a  certain  control  law  and  a  set  of  weighting  matrices.  Whenever 
a  control  mode  is  used  at  different  sampling  rates,  the  gains  are 
recalculated  to  account  for  the  new  sampling  interval.  So,  for  instance. 
Pitch  Rate  Mode  B  for  10  sps  does  not  have  the  same  optimal  gains  as 
Mode  B  for  7  sps  (Table  3-3).  All  plots  are  for  an  initial  condition 
of  straight-and-level  flight  at  105  KIAS.  For  each  run,  the  six  vari¬ 
ables  Aq,  A0,  Anz,  Aa,  A V,  and  A6£  were  plotted.  The  step  input  on 
A6S  was  set  and  then  released  after  approximately  two  seconds. 


> 
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TABLE  3-2 

Control  Modes  Used  in  Hybrid  Tests 


CAS 

Mode 

Q  Weightings 
qll  »  q22 

R  Weightings 

Direct 

. 

Pitch  Rate 

A 

3,20 

33 

11 

B 

25,25 

33 

!» 

C 

50,50 

33 

tt 

D 

75,25 

33 

Norm .  Acc . 

A 

5,10 

33 

It 

B 

10,5 

33 

II 

C 

5,10 

250 

11 
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TABLE  3-3 


Optimal  Gains  for  the  Control  Configurations 


Control 

Mode 

Sampling 

Rate  (sps) 

Gains 

q 

ct(n„) 

q  Mode  A 

10 

-.150 

.036 

q  Mode  B 

10 

-.427 

.324 

q  Mode  C 

10 

-.554 

.360 

q  Mode  D 

10 

-.620 

.469 

q  Mode  C 

8 

-.484 

.377 

q  Mode  D 

8 

-.532 

.475 

q  Mode  B 

7 

-.358 

.347 

q  Mode  B 

6 

-.324 

.355 

q  Mode  B 

5 

-.282 

.362 

q  Mode  B 

4 

-.232 

.366 

n2  Mode  A 

10 

-.533 

(-.181) 

nz  Mode  B 

10 

-.474 

(-.106) 

nz  Mode  C 

10 

-.238 

(-.054) 

nz  Mode  A 

8 

-.494 

(-.154) 
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The  Direct  Mode  (Fig.  3-5a)  demonstrates  the  predicted  open- 
loop  response  of  the  aircraft.  Pitch  rate  overshoots  and  decays, 

Aa  reaches  a  new  steady  state,  and  AV  decreases  as  the  aircraft 
starts  to  climb.  When  Ay<j  is  set  back  to  zero,  Aa  returns  to  its 
original  value; the  decreased  AV  leads  to  a  negative  Aq  and  a  rapidly 
decreasing  A9. 

The  Pitch  Rate  Mode  B  (Fig.  3-Sc)  shapes  the  A<Sg  input,  pro¬ 
viding  quicker  Aq  response  with  less  overshoot.  The  Aq  decay  is  re¬ 
duced  but  not  eliminated  because  the  control  law  is  not  Type  1  and 
there  is  no  throttle  augmentation  to  control  the  decrease  in  AV.  Upon 
releasing  the  stick,  Aq  returns  closer  to  the  original  value  than  the 
previous  case,  and  the  rate  of  A0  decrease  is  reduced.  Comparing 
these  plots  to  those  generated  in  the  analytic  study  (Fig.  2-7),  it 
is  seen  the  Micro-DFCS  is  operating  as  expected.  The  initial  control 
deflection  is  approximately  3.5  degrees  on  both  plots,  and  the  responses 
of  Aq,  Aa,  and  AV  are  similar.  Increasing  the  weighting  on  Q  for  the 
Pitch  Rate  Mode  C  (Fig.  3-5d)  shows  the  expected  slight  decrease  in 
rise  time  for  Aq,  with  about  the  same  amount  of  overshoot.  To  obtain 
the  faster  rise  time,  the  magnitude  of  the  initial  A6E  deflection  has 
increased.  With  Pitch  Rate  Mode  D,  both  rise  time  and  overshoot  have 
decreased  as  Fig.  2-5  predicts.  Different  sampling  rates  were  used 
as  shown  in  Fig.  3-5f-k.  The  optimal  gains  were  recalculated  for 
each  set  of  weightings,  taking  into  account  the  change  in  sampling 
interval.  Comparing  runs  with  the  same  weightings,  but  different 
sampling  rates,  shows  identical  state  trajectories  but  quite  different 
control  inputs,  A6g.  The  state  trajectories  are  the  same  for  different 
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sampling  rates  because,  as  mentioned  in  Chapter  2,  the  cost  functional, 
J,  is  minimized  continuously  in  time  regardless  of  the  sampling  rate. 
The  control  input  is  different  to  account  for  the  different  sampling 
intervals.  In  all  control  configurations,  the  zero-order  hold  steps 
in  AS£  are  visible.  The  state  responses  to  these  control  inputs  are 
smooth  because  the  rigid-body  dynamics  of  the  VRA  filter  out  the 
higher  frequency  content  of  the  sharp-edged  steps. 

The  Normal  Acceleration  Mode  A  (Fig.  3- SI)  forces  a  large 
Aq  overshoot  to  obtain  rapid,  well-damped  An2  response.  The  A6£  de¬ 
flection  begins  with  a  large  impulse;  this  is  followed  by  a  gradual 
deflection  which  increases  Aa  to  minimize  Anz  decay  as  AV  decreases. 
Figure  3-5n  gives  the  response  for  R  =  250.  The  weighting  R  increases 
the  importance  of  suppressing  perturbations  in  Au  in  the  cost  func¬ 
tional  J.  Comparing  Fig.  3-5m  to  Fig.  3-5n,  the  A6£  deflection  is 
initially  less  for  R  =  250,  even  though  the  step  input  used  (Ay,j=.25g) 
is  greater.  Figure  3-5o  was  made  with  the  sampling  rate  set  at  8  sps. 
This  plot  shows  once  again  that  sampling  rate  does  not  effect  the  time 
response  using  sampled-data  regulation  theory. 

The  actual  flight  tests  conducted  for  the  Micro-DFCS  are 
presented  in  the  next  chapter  along  with  the  results  and  analysis. 

This  marks  the  next  major  step  in  the  testing  of  the  implemented 
control  laws  and  the  Micro-DFCS  equipment. 


A«-  (deg)  AV  (fps)  Ad  (deg)  A"z  (g)  A9  (deg)  ^  (de*/s«) 


Figure  3-5. 
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4. 


FLIGHT  TESTING 


Flight  testing  is  an  important  segment  of  the  Micro-DFCS 
implementation.  Experiment  and  demonstration  in  flight  are  necessary 
for  a  full  understanding  of  the  relationships  between  control  theory 
and  practice.  The  qualitative  and  quantitative  results  obtained  in 
this  investigation  show  some  interesting  characteristics  of  DFC  with 
low  sampling  rates  (10  sps  and  belcw).  The  overall  results  indicate 
that  low  sampling  rates  of  the  Micro-DFCS  do  not  cause  destabilizing 
phase  lags  as  would  be  expected  with  digitized  analog  control  laws. 
These  results  are  presented  in  Section  4.2.  Section  4.1  outlines 
the  test  objectives  and  procedures. 

4 . 1  TEST  OBJECTIVES  AND  PROCEDURES 

Two  objectives  were  set  for  the  flight  tests.  The  first 
objective  was  to  check  out  the  Micro-DFCS  installation,  and  its 
operation  with  the  VRA  fly-by-wire  system.  It  was  anticipated  that 
problems  may  have  occurred  due  to 

•  RF  noise 

•  Sensor  and  channel  noise 

•  Scaling  factors  set  incorrectly 

within  the  software 

•  Difficulties  in  using  the  control 

display  unit  while  in  flight 

The  possibility  existed  that  RF  radiation  from  the  telemetry  or  com¬ 
munication  radios  might  effect  the  operation  of  the  Micro-DFCS  or 
vice-versa.  Also  noise  might  be  present  on  the  sensor  signals  (which 
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was  not  accounted  for  in  the  hybrid  simulation)  that  would  degrade 
the  controller's  performance.  Noise  at  high  frequencies  (because  of 
structural  vibrations)  that  would  not  normally  affect  the  performance 
of  the  VRA's  analog  control  system,  could  be  a  problem  with  a  digital 
controller.  The  digital  system  can  cause  the  higher  frequency  con¬ 
tent  of  the  inputs  to  be  "folded"  down  to  the  frequencies  of  interest 
because  of  sampling  effects.  Another  problem  could  be  caused  by  incor¬ 
rectly  set  scaling  factors  that  convert  voltage  levels  of  pilot  con¬ 
trols  and  feedback  sensors  to  the  corresponding  units  of  motion. 
Finally,  there  was  a  question  whether  the  CDU's  display  could  be  seen 
well  enough  in  direct  sunlight,  and  if  the  method  of  keying  in  re¬ 
quests,  defined  in  CAS-1,  was  adequate. 

Examination  of  the  effects  of  changing  gains  and  sampling 
rates  for  the  control  laws  used  in  the  Micro-DFCS  was  the  second 
major  objective  for  testing.  The  influences  of  these  parameters  on 

•  Flying  qualitites  criteria 

•  Pilot's  opinion  of  controllability, 
responsiveness,  and  other  noted  effects 

•  Step  response  as  compared  to  the  ana¬ 
lytic  and  hybrid  results 

•  Steady  turn  maneuvers 

were  specifically  examined.  The  flying  qualitites  criteria  from  the 

Military  Specification  (Ref.  20)  is  given  in  terms  of  and  u>n 

.  nsp 

versus  n2/a.  The  closed-loop  values  of  these  parameters  for  each 
set  of  gains  were  determined  in  the  analytical  study;  the  flight 
tests  were  intended  to  verify  the  flying  qualities  classifica¬ 
tion.  The  pilot's  opinion  of  controllability  8nd  responsiveness 
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is  a  subjective  indication  of  flying  qualities.  The  step  response 
test  was  conducted  with  the  same  step  inputs  as  the  hybrid  test. 

A  step  Ayj  command  was  entered  on  the  CDU  and  the  resulting  motion 
response  was  telemetered  to  the  ground  for  analysis.  The  last  goal 
was  to  insure  that  this  longitudinal  controller  did  not  significantly 
oppose  steady-turn  controller  displacement. 

The  preparations  for  each  flight  test  consisted  mainly  of 
identifying  the  control  configurations  (control  mode  and  sampling 
rate)  to  be  tested  on  each  flight,  and  loading  the  flight  control 
computer  unit  on  its  palet,  which  was  mounted  in  the  VRA  (attached 
with  quick  release  bolts).  The  FCCU  was  reinstalled  for  each  test 
because  CAS-1  usually  needed  slight  modifications  depending  on  the 
control  configuration.  Also  the  Micro-DFCS  equipment  and  software 
were  checked  in  the  hybrid  simulation  prior  to  flight  tests  to  verify 
that  all  scalings  and  gains  were  set  correctly. 

The  testing  procedure  for  each  flight  was  basically  the  same. 
Once  the  VRA  had  attained  the  proper  altitude  (5000  feet  MSL  for  most 
tests),  the  Micro-DFCS  was  engaged  using  the  following  sequence: 

2.  Attain  straight-and-level  flight 
at  105  KIAS 

2.  Select  control  configuration 

3.  Radio  type  of  test  to  ground  station 

4.  Engage  VRA’s  FBW  system 

5.  Accomplish  test  maneuver 

6.  Disengage  FBW  system 

7.  Set  up  for  next  test 

The  first  step  is  necessary  to  put  the  VRA  in  the  control  design  flight 
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condition.  Some  tests,  however,  were  designed  to  find  how  the  con¬ 
troller  operates  when  under  different  flight  conditions;  in  that 
case  Step  1  was  modified  accordingly.  The  next  step  selects  the 
control  law,  gains,  and  sampling  rate  by  entering  the  requests  through 
the  CDU  (see  Appendix  D  for  a  detailed  description  for  using  the  CDU). 
When  a  new  control  mode  (different  control  law  or  set  of  gains)  is 
requested  using  CAS-1,  the  feedback  states  at  that  instant  are  taken 
as  the  "set  point"  or  nominal  value  about  which  the  Micro-DFCS  con¬ 
trols  for  zero  A6g.  Therefore,  it  is  important  that  the  flight  con¬ 
dition  desired  for  the  run  be  maintained  while  Step  2  is  being  performed. 
Step  3  lets  the  ground  station  operator  mark  the  strip  chart  plots  and 
tape  recording  with  the  type  of  run.  It  is  done  just  prior  to  the 
maneuver  so  the  strip  chart  and  tape  recorder  do  not  run  needlessly 
during  the  set-up.  Next,  the  FBW  system  is  engaged,  allowing  the  6E 
output  from  the  Micro-DFCS  to  operate  the  elevator  actuator.  Up  to 
this  point  in  the  sequence,  the  Micro-DFCS  calculates  a  <5^  but  the 
actuator  is  not  enabled  to  follow  the  signal.  When  the  FBW  system 
is  engaged,  the  Micro-DFCS  6g  signal  is  automatically  biased  so  that 
it  corresponds  to  the  actual  elevator  position  just  before  the  engage¬ 
ment.  This  insures  a  smooth  transition  from  the  normal  mechanical 
controls  to  the  FBW  controls.  The  maneuver  is  performed  and  the  test 
data  for  that  run  then  is  collected.  The  FBW  system  is  disengaged, 
and  the  sequence  is  repeated  for  the  next  test  run. 

The  tests  yielded  both  qualitative  and  quantitative  results. 

In  one  set  of  tests,  the  pilot  performed  a  defined  task  and  rated 
the  responsiveness  and  controllability  for  each  control  configuration. 
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Quantitative  tests  were  conducted,  where  the  response  of  the  VRA  to 
a  step  input  on  <$s  was  telemetered  for  analysis.  Figure  4-1  shows 
the  breakdown  of  the  test  methods  used. 


Figure  4-1.  Breakdown  of  Test  Methods  Used. 

The  tracking  task  was  conducted  by  engaging  the  Micro-DFCS  and  then 
attempting  to  position  a  "pip"  or  mark  on  the  windshield  over  an 
arbitrary  point  on  the  horizon.  Every  few  seconds  another  point  on 
the  horizon  was  picked  and  the  pilot’s  success  in  tracking  the  new 
point  was  recorded  as  the  results  for  that  test.  The  pilot  also  was 
asked  to  comment  on  his  impression  of  the  discrete  elevator  movements. 
The  tracking  task  was  performed  for  the  Direct  Mode  and  Pitch  Rate 
Mode  B  at  sampling  rates  of  10,  8,  7,  6,  5,  and  4  sps.  It  was  also 
done  for  all  Normal  Acceleration  Modes  and  Pitch  Rate  Modes  C  and  D 
at  10  sps. 
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The  landing  approach  test  was  conducted  only  in  the  Direct 
Mode  because  both  control  inodes  were  developed  only  for  a  105  KIAS 
flight  condition.  Only  the  effects  of  the  sampling  delay  were  being 
examined  in  the  landing  test.  The  testing  period  for  each  landing 
started  right  after  the  turn  to  the  final  approach  until  touchdown. 
Sampling  rates  of  10,  7,  5,  3,  2,  and  1  sps  were  used. 

Analysis  of  the  response  of  the  commanded  state  (either  Aq 
or  Anz)  to  a  step  input  of  A6g  was  examined  as  in  the  hybrid  tests. 

The  same  step  inputs  were  used  so  a  comparison  could  be  made.  All 
control  configurations  presented  in  Section  3.3  were  tested  in  flight. 
In  addition,  the  Pitch  Rate  Mode  B  operating  at  10  sps  was  tested  in 
straight-and- level  flight  at  80  KIAS  and  again  at  120  KIAS. 

The  last  testing  method  was  performed  to  analyze  the  effects 
of  this  longitudinal  controller  when  attempting  to  hold  a  steady  co¬ 
ordinated  turn  with  bank  attitudes  (<j>)  of  30  and  50  deg.  The  test 
was  done  using  the  Direct  Mode  and  the  Pitch  Rate  Mode  B  at  10  sps. 

The  questions  to  be  answered  for  this  test  were: 

1.  How  large  an  effect  does  the  pitch  rate 
damping  have  on  steady-state  control  position? 

2.  Is  this  deflection  objectionable  to  the  pilot? 

4.2  FLIGHT  TEST  RESULTS  AND  ANALYSIS 

Three  flights  were  made  to  collect  the  data.  The  first  flight 
was  used  to  test  the  components  of  the  VRA/Micro-DFCS  combination. 

The  components  did  not  interfere  with  each  other's  operation  except  in 
one  case.  Whenever  the  telemetry  transmitter  was  turned  off,  the  FCCU 


4-6 


4 


of  the  Micro-DFCS  was  halted.  It  was  not  clear  if  this  was  due  to 
an  RF  pulse  or  a  pulse  on  the  power  lines.  This  did  not  create  a 
problem  for  testing  the  control  laws,  so  it  was  not  corrected. 

The  general  operation  of  all  but  the  Normal  Acceleration  Mode 
was  tested  on  the  first  flight  because  the  nz  accelerometer  signal  was 
too  noisy  to  allow  the  control  law  to  function.  Stick-gearing  for 
the  Pitch  Rate  Modes  was  found  to  be  too  low  for  effective  handling 
qualities  evaluation.  It  was  originally  set  to  command  2.7  deg/sec 
per  inch  of  stick  deflection  in  the  steady-state.  The  pilot  considered 
this  too  much  movement  of  the  stick  for  the  desired  Aq.  The  stick¬ 
gearing  was  doubled  which  gave  adequate  control  power  at  the  105  KIAS 
flight  condition.  The  Normal  Acceleration  Mode  stick-gearing  was  also 
doubled  from  its  previous  value  to  0.52  g/inch.  The  stick-gearings 
were  set  low  on  Flight  1  for  overflow  protection;  for  some  control  modes 
a  large  stick  deflection  can  cause  the  Micro-DFCS  to  temporarily  com¬ 
mand  more  than  the  maximum  6g  allowable.  This  causes  the  D/A  output 
to  overflow,  resulting  in  a  jump  in  output  voltage,  as  described  for 
the  LIMIT  routine  in  Appendix  E.  The  routine  LIMIT  was  added  to  CAS-1 
after  Flight  1  to  prohibit  D/A  overflow. 

Noise  on  the  sensor  channels  was  a  problem  for  the  nz  accel¬ 
erometer.  Structural  vibrations  with  frequencies  of  10-Hz  or  more 
were  superimposed  on  the  nz  signal.  A  first-order  low-pass  analog 
filter  with  4-Hz  bandwidth  was  added  to  the  n?  channel  before  the  A/D 
converter  input.  For  the  short-period  natural  frequency  (fSp«0. 56Hz) , 
this  filter  causes  a  one  percent  attenuation  and  an  8  deg  phase  lag 
in  the  nz  signal.  The  filter  completely  eliminated  the  effects  of 


the  structural  vibration  on  the  Micro-DFCS  in  the  remaining  test 
flights. 
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The  second  and  third  flights  tested  the  modifications  made 
after  Flight  1  and  provided  data  for  the  qualitative  and  quantitative 
tests  outlined  in  Section  4.1. 

Tracking  at  Altitude 

The  results  for  the  tracking  task  are  given  in  Table  4-1. 

The  pilot's  comments  for  each  control  mode  and  sampling  rate  used  are 
listed.  The  temperature  was  17°C,  the  altitude  was  4500  to  5500  ft  MSL, 
and  barometric  pressure  was  30.16  inches.  The  tests  were  conducted  in 
essentially  calm  air.  Normal  Acceleration  Mode  A  was  expected  to  give 
an  undesirable  response  because  the  closed- loop  for  this  mode 

(^^=8.85  rad/sec,  5Sp=.73)  places  the  flying  qualities  in  the  Level 
2  region  of  Fig.  2-3.  Figure  2-3  is  presented  again  in  Fig.  4-2  with 
all  control  modes  marked.  Normal  Acceleration  Mode  B  is  also  in  Level 
2  with  (%  =  7.62  rad/sec.  The  flying  qualities  on  n  Mode  C  is  in 

Level  1  with  =  5.5  rad/sec;  this  rating  is  also  reflected  in  the 

pilot's  comments.  The  discrete  pulsing  movements  of  the  elevator 
can  be  felt  and  heard  in  varying  degrees,  depending  on  the  sampling 
rate  and  control  mode.  Low  sampling  rates  or  control  modes  that  com¬ 
mand  large  initial  deflections  make  the  discrete  nature  of  the  move¬ 
ments  very  noticeable.  The  vibrations  seem  very  unnatural  and  are 
likened  to  the  sound  and  feel  of  a  stall  buffet.  The  phenomena  is 
very  annoying  and  a  DFCS  with  zero-order  hold  outputs  at  less  than 
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Figure  4-2.  Short-Period  Frequency  Requirements 
with  All  Control  Modes  Marked. 


10  sps  could  probably  not  be  tolerated  in  an  operational  design  for 
this  flight  condition  using  the  control  structure  defined  in  Section 
2.1;  however,  no  conclusion  can  be  reached  regarding  the  suitability 
of  these  low  sampling  rates  with  more  advanced  control  designs  that 
smooth  control  rates. 


The  Direct  Mode  operating  at  10  sps  gives  the  normal  open- 
loop  response  of  the  aircraft  with  no  noticeable  sampling  effects. 


TABLE  4-1 


Results  for  Tracking  Task 


Sampling 
Rate  (sps) 


Control 

Mode 

nz  Mode  A 


nz  Mode  B 


n2  Mode  C 


Direct  Mode 


q  Mode  B 


q  Mode  C 
q  Mode  D 


Direct  Mode 


q  Mode  B 


Pilot's 

_ Comments _ 

Movements  too  abrupt.  Over  sensi¬ 
tive  to  command  inputs.  Large 
initial  6g  pulses  are  felt  and 
heard.  FBW  system  disengages  on 
large  input  commands. 

Better  than  Mode  A,  but  initial 
response  still  toe  fast.  Discrete 
deflections  still  noticeable. 

Good  response  witl  "crisp"  track¬ 
ing.  Discrete  movements  sound 
and  feel  like  the  stall  buffet. 

This  seems  very  unnatural. 

No  noticeable  difference  from  the 
normal  continuous  open- loop 
response  except  for  slight  rum¬ 
ble  of  6g  deflection. 

Good  response.  Easier  to  attain 
track  of  mark  and  hold  then  Direct 
Mode.  "Rumble"  from  elevator 
about  same  as  Direct  Mode. 

Seems  the  same  as  q  Mode  B. 

Response  and  control  the  same  as 
q  Modes  B  and  C,  but  rumble  from 
elevator  more  noticeable.  Not 
really  distracting,  but  obviously 
present . 

The  same  response  as  Direct  Mode 
for  10  sps.  Rumble  may  be  slightly 
more  apparent . 

The  same  response  as  q  Mode  B  for 
10  sps.  Rumble  may  be  slightly 
more  apparent. 


(continued) 


TABLE  4-1. 


(continued) 


Sampling 
Rate  (sps) 

Control 

M  ode 

Pilot's 

Comments 

7 

Direct  Mode 

A  little  overshoot  noticed  on 
tracking  task  but  does  not  inter¬ 
fere  with  performance  significantly. 
Rumble  of  6g  movements  more  notice¬ 
able. 

7 

q  Mode  B 

No  overshooting  of  mark.  Rumble 
as  in  the  Direct  Mode  above. 

6 

Direct  Mode 

Continual  overshoot  of  mark  in 
tracking  task.  Very  objectionable, 
but  tolerable. 

6 

q  Mode  B 

Overshoot  of  mark  noticed,  but 
tracking  attained  more  quickly  than 
Direct  Mode  at  6  sps.  Rumble  now 
more  of  a  "thumping";  more  notice¬ 
able  than  Direct  Mode. 

5 

Direct  Mode 

Response  seems  better  than  Direct 
or  q  Mode  at  6  sps.  Thumping 
very  noticeable  (sound  and  vibra¬ 
tions. 

5 

q  Mode  B 

Worse  than  Direct  Mode.  Tendency 
to  overshoot  mark,  but  can  even¬ 
tually  get  there.  6g  thumping 
about  the  same  as  above. 

4 

Direct  Mode 

Performance  very  bad.  Very 
large  overshoots  in  tracking 
task. 

4 

q  Mode  B 

Pilot  senses  delay  from  command 
until  motion,  but  no  trouble  in 
attaining  and  keeping  the  track. 

No  overshoot,  just  a  very  objec¬ 
tionable  delay.  thumping  same 

as  5  sps. 
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The  elevator  movements  are  noted  as  a  kind  of  "rumbling"  sound.  The 
Pitch  Rate  Modes  B  and  C  at  10  sps  give  a  better  response  than  the  Direct 
Mode  because  the  overshoots  and  rise  times  in  q,  arising  from  short-period 
modal  characteristics,  are  decreased.  Pitch  Rate  Mode  D  gives  a  similar 
response,  but  initial  elevator  deflections  are  larger  with  Mode  D, 
which  causes  the  deflections  to  be  more  noticeable.  The  flying  qualities 
criterion  set  by  Fig.  4-2  shows  Level  2  response  for  Modes  C  and  D  with 

=  6.5  rad/sec,  and  7.2  rad/sec  respectively  (taken  from  Table  2-2). 
sp 

The  pilot's  comments  indicate,  however,  that  the  response  and  control¬ 
lability  are  very  good  for  these  two  modes.  An  explanation  may  be 
that  when  the  £Sp  is  very  high  (in  this  case  1.5  and  1.8),  Fig.  4-2 
may  not  necessarily  apply.  The  Military  Specification  gives  no  neces¬ 
sary  relationship  between  uu  and  t,  ;  possibly  one  is  needed. 

sp  bP 

The  remaining  tests  listed  in  Table  4-1  are  for  sampling 
rates  8  through  4  sps,  comparing  the  response  for  the  Direct  Mode 
and  Pitch  Rate  Mode  B.  As  the  sampling  rate  is  decreased,  the  ele¬ 
vator  vibrations  get  worse  and  are  described  as  a  "thumping"  sound. 

The  deflections  can  also  be  felt.  The  Pitch  Rate  Mode  is  always 
worse  in  this  respect  than  the  Direct  Mode  (at  the  same  sampling  rate). 
Down  to  8  sps,  the  response  is  considered  the  same  as  for  10  sps  for 
both  the  Direct  Mode  and  Pitch  Rate  Mode  B.  At  7  sps  and  below,  a 
definite  degradation  in  performance  is  noted  in  the  Direct  Mode  because 
of  a  tendency  for  the  pilot  to  overshoot  the  mark.  This  overshoot 
probably  is  caused  by  the  pilot  expecting  an  immediate  pitching  motion 
that  is  instead  delayed  by  the  sampling  interval.  When  the  anticipated 
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motion  does  not  occur,  the  pilot  moves  the  stick  slightly  more,  which 
causes  too  much  when  the  sampling  instant  does  occur.  Each  Pitch 
Rate  Mode  B  is  better  than  the  Direct  Mode  for  the  same  sampling 
rate,  except  for  the  5  sps  case,  a  seemingly  anomalous  result.  At 
4  sps,  the  pilot  definitely  can  sense  a  delay  but  comments  that  he 
has  no  trouble  attaining  and  holding  a  track  on  the  mark.  This  could 
be  a  learning  effect  which  occurred  during  these  particular  tests. 

The  pilot  may  control  overshoot  better  with  Pitch  Rate  Mode  B  at 

4  sps  because  he  finally  realizes  why  he  is  overshooting  the  mark 
when  a  definite  delay  is  sensed.  The  pilot  may  not  recognize  the 
problem  with  the  Direct  Mode,  because  the  overshoots  in  Aq  could  be 
attributed  to  the  lighter  damping  of  the  unaugmented  short  period  mode. 

As  mentioned  above,  the  Direct  Mode,  for  some  reason,  seems 
better  than  the  Pitch  Rate  Mode  at  the  5  sps  rate.  The  5  sps  tests 
were  repeated  to  assure  that  no  error  in  set-up  was  made,  each  time 
with  the  same  result.  No  explanation  has  been  found  for  the  degraded 

5  sps  results.  Possible  (unverified)  causes  include  calculation  of 
gains  done  incorrectly,  or  that  there  is  some  type  of  interaction 
between  the  sampling  frequency  and  the  frequency  of  the  pilot's  inputs. 

Landing  Approach 

The  results  for  the  landing  approach  test  are  given  in  Table 
4-2.  The  atmospheric  conditions  are  the  same  as  for  the  tests  listed  in 
Table  4-1.  The  tests  were  conducted  in  essentially  calm  air.  Analysis 
of  this  test  shows  the  flare  maneuver  to  be  less  sensitive  to  sampling 
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TABLE  4-2 


Results  for  Landing  Approach  Using  the  Direct  Mode 


Sampling 

Rate  (sps) 

Pilot's  Comments 

10 

No  problem:  not  aware  of  sampling  on 
vibrations. 

7 

Slight  overshoot  in  response  on  short 
final,  but  no  problem  in  flare. 

5 

Larger  overshoots  on  short  final,  but 
still  no  problem  in  flare. 

3 

Overshoots  very  objectionable  on  short 
final,  flare  still  no  problem. 

2 

The  marginal  case:  overshoots  cause 
poor  control  on  both  short  final  and 
flare. 

1 

Not  possible  in  this  configuration: 
pilot  moves  stick  far  enough  in  1  sec 
to  send  a  step  to  that  disengages 

the  safety  feature  of  the  FBW  system  of 
VRA. 

rate  than  the  short  final  approach  or  tracking  task.  This  might  be 
explained  by  examining  the  pilot's  reaction  to  delayed  response.  In 
the  tracking  task  at  10S  KIAS  or  on  short  final  at  70-75  KIAS,  the 
control  surfaces  generate  a  larger  moment  than  in  the  slower  flare 
(60  KIAS).  The  pilot,  therefore,  expects  an  immediate  response  to 
his  commanded  inputs  when  at  cruising  speed  or  on  the  short  final 
approach,  but  he  has  learned  there  will  be  slower,  more  sluggish 
movements  in  the  flare.  He  anticipates  a  natural  delay  in  response 
on  the  flare  which  in  effect  "masks"  the  delay  of  the  sampling  rate. 


When  the  aircraft  does  not  move  immediately  in  response  to  a  commanded 
input,  the  pilot  expects  this  and  does  not  over  correct  with  a  larger 
command . 


Pitch  Rate  in  Steady  Turn 

The  telemetry  records  for  the  banked  turn  test  shows  that 
the  longitudinal  controller  does  oppose  the  execution  of  a  coordinated, 
banked  turn.  The  results  are  given  in  Table  4-3. 


TABLE  4-3 


Results  for  Banked  Turn  Test 


Control 

Mode 

Bank  Angle 
(deg) 

Elevator 

Deflection 

(deg) 

Steady  Longitudinal  Stick 
Deflection  Required,  (inches) 

Direct 

30 

-1 

.24 

q  Mode  B 

30 

-1 

.5 

Direct 

50 

-3 

1.2 

q  Mode  B 

50 

-3 

1.7 

More  6S  is  needed  for  the  Pitch  Rate  Mode  than  the  Direct  Mode  at 
both  bank  angles.  This  is  true  because  the  Micro-DFCS  pitch  rate 
feedback  is  trying  to  suppress  the  positive  Aq  developed  in  a  steady 
turn.  With  the  stick-force  per  g  characteristic  used  in  these  tests, 
the  pilot  was  not  aware  of  the  additional  deflection  needed,  so  the 
longitudinal  controller  did  not  seem  to  degrade  lateral-directional 
maneuvering  of  the  VRA  in  any  significant  way. 
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The  final  results  obtained  from  the  flight  tests  are  the 
responses  to  commanded  step  inputs.  The  telemetry  records  of  the 
step  responses  are  presented  in  Fig.  4-3.  Because  of  a  technical  mal¬ 
function  of  the  ground  station,  many  of  the  step  response  runs  were 
not  recorded  on  tape  correctly  to  be  played  back  for  better  analysis. 
The  only  telemetry  records  available  are  those  that  came  off  the 
strip  chart  recorder  at  the  time  of  the  actual  test.  The  sense  of 
these  plots  is  opposite  to  that  of  the  analytical  and  hybrid  tests. 

The  axes  are  labeled  correctly,  so  comparisons  still  can  be  made. 

The  analysis  of  the  step  responses  from  the  flight  and 
hybrid  tests  is  done  in  the  same  manner.  Since  there  is  a  corre¬ 
sponding  hybrid  plot  for  each  flight  test  plot  (except  for  the  120 
and  80  KIAS  cases),  and  the  hybrid  plots  were  analyzed  in  Section  3.3, 
only  the  different  or  important  aspects  of  the  flight  test  plots  are 
examined  here.  Pitch  Rate  Modes  B,  C,  and  D  (Fig.  4-3c,d,e)  at  10  sps 
improve  the  Aq  response  essentially  the  same  amount.  Mode  D  commands 
a  larger  initial  A6£  deflection,  as  expected.  As  the  sampling  rate 
is  decreased  for  Pitch  Rate  Mode  B  (Fig.  4-3g,h,i, j ,k) ,  the  Aq  re¬ 
sponse  stays  the  same,  but  the  initial  ASg  deflections  are  smaller, 
to  account  for  the  decrease  in  sampling  rate.  In  Fig.  4-3k  it  is 
easy  to  see  the  0.25  sec  sampling  interval  in  A6E.  Comparing  plots 
from  Fig.  4-31  and  m  to  Fig.  4-3c,  there  is  not  much  difference  in  Aq 
response  from  speeds  of  80  to  120  KIAS.  At  80  KIAS  it  is  apparent, 
however,  that  more  A6^  is  needed  to  hold  the  commanded  Aq  as  velocity 
falls  to  even  lower  values. 


An  important  difference  between  the  Pitch  Rate  Mode  plots 
of  the  hybrid  and  flight  tests  is  the  presence  of  small  disturbances 
on  the  latter.  It  is  probably  not  noise  generated  in  the  FBW  or 
telemetry  systems  because  the  Normal  Acceleration  Mode  responses 
are  relatively  smooth  (Fig.  4-3n) .  Most  likely  it  is  small  air  dis¬ 
turbances  at  the  wing  tips  picked  up  by  the  a  vanes.  These  dis¬ 
turbances  can  be  seen  in  the  Aa  plots  in  Fig.  4-3g  and  m.  The  a 
signal  feeds  to  the  elevator,  which  in  turn  causes  the  small  pertur¬ 
bations  in  Aq.  This  is  a  good  example  of  why  air  data  sensors  might 
not  be  the  best  choice  for  some  designs.  In  this  case,  however,  the 
pilot  is  not  aware  of  these  small  perturbations  and  gives  the  Pitch 
Rate  Modes  B-D  good  ratings  at  10  sps. 

The  Normal  Acceleration  Modes  have  plots  in  Fig.  4-3  similar 
to  the  hybrid  plots.  Normal  Acceleration  Mode  A  has  an  initial  A6g 
deflection  that  is  too  large  to  use  a  commanded  input  of  yj  *  .25g. 

To  test  Normal  Acceleration  Mode  A,  a  commanded  input  of  y^  *  .17g 
is  used  to  keep  the  FBW  system  from  disengaging.  The  large  overshoots 
in  Aq  are  present  as  expected. 

The  next  and  final  chapter  summarizes  the  major  results  of 
all  the  previous  chapters.  Recommendations  are  made  for  continuation 
of  research  in  this  project. 
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Figure  4-3.  continued 
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CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

The  results  of  this  study  provide  a  basis  for  low-cost 
flight  testing  of  advanced  concepts  in  digital  flight  control  using 
both  modern  control  theory  and  microprocessor  technology.  The 
proof-of-concept  control  laws,  derived  using  optimal  control  theory, 
allow  very  low  sampling  rates  to  be  used  (10  sps  or  less).  This  is 
because  the  gains  for  the  simple  Type  0  pitch  rate  and  normal  ac¬ 
celeration  CAS  structures  are  derived  by  minimizing  an  analog  per¬ 
formance  index  even  though  the  control  law  is  calculated  at  discrete 
instances  in  time.  This  gives  a  large  improvement  in  sampling  re¬ 
quirements  over  the  conventional  method  of  simply  digitizing  an  ana¬ 
log  control  law.  The  resulting  lower  sampling  rates  decrease  the 
computational  burden  associated  with  the  controller  and  would  allow 
more  advanced  control  modes  to  be  implemented  on  the  microprocessor 
based  system.  The  Micro-DFCS  is  integrated  with  the  VRA;  the  system 
is  completely  operational  and  ready  for  further  digital  flight  con¬ 
trol  research.  New  flight  test  programs  can  be  developed  quickly 
because  all  hardware  is  installed,  and  a  basis  for  further  software 
development  is  present  in  CAS-1. 

The  objectives  of  this  project,  as  outlined  in  Chapter  1 
were  to  add  a  Micro-DFCS  to  the  VRA's  FBW  analog  control  system 
and  provide  initial  results  with  actual  flight  testing.  The 


objectives  have  been  met  with  all  results  and  analysis  given  in  the 
text  of  this  report.  A  summary  of  the  major  conclusions  follows. 

•  The  step  response  plots  from  the  digital  model, 
analog  simulation,  and  flight  tests  all  show  that  the 
major  response  features  of  the  closed-loop  system  are 
not  affected  by  sampling  rate  when  sampled-data  regu¬ 
lator  theory  is  used.  This  verifies  that  the  lag 
which  results  from  sampling  the  feedback  variables 

is  indeed  accounted  for  in  the  control  law  calcula¬ 
tion.  The  lower  limit  on  sampling  rate  in  this 
study  is,  therefore,  not  determined  by  the  phase 
lag  inherent  in  the  sampling  process. 

•  The  microcomputer  used  in  this  study  has  the 
speed  to  execute  much  more  complicated  and  advanced 
control  laws  at  the  sampling  rates  examined.  For 
example.  Type  0  fourth-order  longitudinal  and 
fourth-order  lateral-directional  control  laws  re¬ 
quire  84  calculations  per  sampling  interval  com¬ 
pared  to  20  for  the  two-state  feedback  control 
laws  examined  here  (including  all  scalings). 

This  implies  that  the  control  laws  for  an  eight-state 
feedback  longitudinal,  lateral-directional  control¬ 
ler  can  be  calculated  in  less  than  26  msec  (84/20 
x  6  msec).  For  a  nominal  sampling  rate  of  10  sps, 
over  74%  of  the  available  duty  cycle  still  is  not 
being  used. 
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•  Commercially  available  microcomputer  equipment 
provides  substantial  capability  for  conducting  ad¬ 
vanced  research  in  DFC.  Additional  memory,  analog 
channels,  or  other  peripheral  devices  can  be  added 
quickly  to  the  system.  The  ease  of  this  expansion, 
coupled  with  the  speed  of  the  machine,  makes  it 
very  versatile  for  further  research. 

•  The  software  development  system  provides  an 
efficient  means  of  writing  and  testing  the  Micro- 
DFCS  software.  The  rudimentary  text  editor  of 
the  microcomputer  allows  easy  entry  of  assembly 
language  programs.  Editing  the  programs  on  the 
microcomputer  is  a  cost-effective  alternative  to 
using  the  time-sharing  computer  for  this  purpose. 

•  Modifying  the  commanded  state's  response  by 
changing  the  weighting  factors  in  the  Q  and  R 
matrices  is  not  difficult,  but  it  is  an  iterative 
procedure,  and  it  is  not  always  obvious  how  the 
weightings  should  be  changed  to  attain  the  de¬ 
sired  response.  Figures  2-S  and  2-12  show  the 
relationships  between  weightings  and  response  for 
the  particular  CAS  structures  of  this  study,  but, 
in  general,  there  is  no  standard  relation  between 
step  response  and  quadratic  weighting  factors. 

These  figures  may,  however,  give  some  indication 
of  what  can  be  expected  for  similar  control  law 
developments. 


•  The  pilot  is  not  aware  of  the  sampled  nature  of 
the  Micro-DFCS  at  10  sps  for  the  Pitch  Rate  Modes 
and  Direct  Mode.  At  this  rate,  no  effects  on  flying 

qualities,  due  to  sampling,  are  seen,  and  the  sound 
and  vibrations  from  the  step  elevator  movements  are 
barely  noticeable.  For  the  pitch  rate  CAS  at  105 
KIAS ,  performance  degradation  is  not  noted  down  to 
7  sps;  however,  the  elevator  movements  are  rather 
annoying.  They  sound  and  feel  similar  to  the  cues 
sensed  in  a  stall  buffet.  A  sampling  rate  of  3  sps 
is  found  to  be  acceptable  using  the  Direct  Mode  on 
the  flare  of  the  landing,  but  a  minimum  of  10  sps 
is  needed  for  the  short  final  segment  of  the  approach. 

•  The  nz  response  cannot  be  speeded  up  effectively 
at  the  105  KIAS  flight  condition  with  only  an  ele¬ 
vator  control  input.  For  any  appreciable  decrease 
in  nz  rise  time,  the  closed-loop  <%Sp  increases  to 
a  value  which  places  the  response  of  the  system  in 
the  Level  2  category.  The  pilot  perceives  the  re¬ 
sponse  as  being  too  touchy  and  oversensitive.  If 
n2  response  is  to  be  speeded  up  at  this  flight 
condition,  direct  lift  control  is  needed. 

The  Micro-DFCS,  using  modern  control  theory,  definitely  has 
a  place  in  future  flight  control  applications  because  of  its  low-cost 
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and  flexibility.  The  modern  control  theory  allows  more  computa¬ 
tions  per  sampling  interval  which,  in  turn,  implies  a  capability  to 
execute  complex  control  functions.  This  increases  the  desirability 
of  the  DFCS  and  especially  a  low-cost  Micro-DFCS.  In  this  way, 
microprocessor  technology  and  modern  control  theory  will  comple¬ 
ment  one  another  and  promote  each  other's  use. 

5.2  RECOMMEND AT IONS 

From  lessons  learned  in  this  study,  five  recommendations 
are  made  for  continuing  research  in  microprocessor  based  digital 
flight  control: 

•  At  sampling  rates  below  8  sps,  some  type  of  out¬ 
put  to  the  control  effectors,  other  than  zero-order 
hold,  may  be  needed  if  control  laws  with  control 
rate  restraint  are  not  used.  The  sound  and  the  vi¬ 
brations  are  too  objectionable  at  the  lower  sampling 
rates.  Possible  corrections  of  the  problem  could 
include  low-pass  filtering  of  the  controller  output 
or  using  output  rates  higher  than  the  input  rates. 

It  must,  however,  be  remembered  that  the  sampled- 
data  regulator  theory  used  in  this  study  assumes 
zero-order  hold  on  the  controller  output  when  mini¬ 
mizing  the  analog  performance  index,  J. 
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•  It  may  be  necessary  to  filter  the  signals  from 
the  sensors  or  pilot’s  controls  before  the  A/D 
conversion  to  prevent  higher  frequency  noise 
from  effecting  the  performance  of  the  Micro-DFCS. 

The  sampling  process  folds  the  higher  frequency 
content  of  the  signal  onto  the  lower  frequencies 

of  interest  and  can  make  the  sampled  signal  useless 

if  the  higher  frequency  content  is  large  in  magnitude 
compared  to  the  low  frequency  information.  For  an 

operational  design,  all  controller  inputs  would 
probably  need  to  be  filtered  to  guard  against 
noise  that  may  develop.  Only  the  n2  accelerometer 
was  filtered  in  this  investigation. 

•  Some  provision  for  trimming  would  be  very  use¬ 
ful  and  improve  efficiency  in  flight  tests.  Once 
the  fly-by-wire  system  is  engaged  in  the  VRA,  there 
is  presently  no  digital  mechanism  to  change  the  set 
point  about  which  the  Micro-DFCS  regulates  with  the 
simple  system  of  this  study.  This  means  either 
that  the  FBW  system  must  be  disengaged  and  controls 
reset  everytime  the  aircraft  drifts  from  the  desired 
flight  condition  or  that  the  VRA' s  analog  system 
must  be  used  to  achieve  trim.  The  slow  integration 
of  a  trim  switch  deflection,  summed  into  y^  may 
provide  the  necessary  control. 
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•  The  Flight  Control  routines  should  be  written 

to  conserve  time  and  not  memory.  The  computational 
speed  is  very  important,  so  the  least  number  of 
program  instructions  should  be  executed  on  each 
calculation.  This  implies  avoiding  the  use  of 
generalized  sub-routines,  and  making  the  Flight 
Control  routines  very  specific  to  their  function. 

•  The  battery  back-up  memory  board  has  been  a  big 
a?*;et  in  the  development  and  use  of  the  microcomputer 
software.  The  transfer  of  programs  from  the 
software  development  system  to  the  FCCU  is 
accomplished  very  easily,  and  program  modifi¬ 
cations  are  much  simpler  and  faster  than  would 

be  the  case  using  PROM.  For  routines  that  sel¬ 
dom  need  to  be  changed,  PROM  storage  would  be 
the  cheapest  and  most  convenient  method.  For 
this  initial  programming  effort,  more  battery 
back-up  memory  will  benefit  the  research  pro¬ 
gram. 

There  are  many  other  variations  and  capabilities  of  the 
Micro-DFCS  that  need  to  be  investigated.  A  few  that  were  recognized 
from  the  results  of  this  study  are  given  below. 

•  A  lower  limit  on  sampling  rate  needs  to  be  esta¬ 
blished  in  the  presence  of  worst-case  turbulence. 
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The  limits  on  sampling  rate  given  in  this  study 
are  set  only  by  the  pilot's  acceptance  criteria 
explained  in  Section  2.2.  The  limit  set  by  error 
build-up  remains  to  be  found. 

•  The  pilot's  comments  on  Pitch  Rate  Mode  D 
show  it  to  exhibit  good  response,  yet  the  Military 
Specification  on  Flying  Qualitites  (Ref.  19)  rates 
this  modi;  in  the  Level  2  category  because  of  the  high 
Wpsp.  This  mode  has  a  rather  large  £sp  which  might 
indicate  that  the  military  specification's  criterion  is 
not  true  at  high  values  of  £sp.  Some  relation 

may  be  needed  in  the  specification  between  u^p 
and  ^sp  ma^e  the  criterion  valid. 

•  Because  of  the  zero-order  hold  on  the  controller 
outputs,  the  control  actuators  always  are  commanded 
to  move  in  steps,  at  infinite  rates.  The  question 
is  raised  as  to  what  kind  of  stress  this  places  on 
the  actuator  mechanism  when  it  always  is  moved  at 
its  maximum  rate. 

•  In  the  tracking  task  of  the  flight  tests,  the 
pilot  preferred  5  sps  in  the  Direct  Mode  to  5  sps 
in  the  Pitch  Rate  Mode.  This  is  contrary  to  the 
preference  at  other  sampling  rates.  While  it  is 


possible  that  a  mistake  was  made  in  calculating  or 
entering  the  gains  in  CAS-1,  the  step  response 
test  for  the  Pitch  Rate  Mode  at  5  sps  was  accom¬ 
plished  during  the  same  flight,  with  the  same 
gains,  and  those  results  are  as  expected  (Fig.4-3j). 
Further  investigation  of  this  anomaly  is  warranted. 


APPENDIX  A 


DERIVATION  OF  DISCRETE  WEIGHTING  MATRICES  FROM  CONTINUOUS 
IMPLICIT  MODEL-FOLLOWING  WEIGHTING  MATRICES 


The  discrete  weighting  matrices,  Q,  M,  and  R  derived  in 
Ref.  15  use  the  continuous  cost  functional  of  Eq.  23.  The  implicit 
model-following  presented  in  Section  2.6  calculates  Q,  M,  and  R 
matrices  for  the  continuous  cost  functional  of  Eq.  56.  In  order 
to  use  the  implicit  model-following  technique  for  a  digital  system, 

A  A  A 

a  transformation  from  Q,  M,  and  R  to  Q,  M,  and  R  is  needed.  This 
relation  is  derived  below  with  the  same  assumptions  used  for  the 
derivations  of  Ref.  12. 

The  objective  is  to  minimize  the  continuous-time  cost  func 
tional,  so  it  is  necessary  to  find  Q,  M,  and  R  that  satisfy  the 
relation 


[xT(t)  Q  x(t)  ♦  2  uT(t)  M  x(t)  ♦  uT(t)  R  u(t) ]dt 


V(*1TQ*i.2ulMxi.ujR»i) 


N  =  number  of  samples 
i  =  index 


given  the  system  dynamics,  Q,  M,  and  R.  The  continuous  J  can  be  de 
composed  into  a  sum  of  N  integrals: 


N-i  1  T  T  T 

J  *  E  (x  Q  X  +  2  u  M  x  ♦  u  R  u)  dt 

i*0  J 


where  the  tine  dependency  is  understood.  Setting  Eq.  64  equal  to  the 
bottom  line  of  Eq.  63  and  making  a  substitution  with 


x(t)  =  flKt.tj)  +  Tft.tj)  Uj 


the  derivation  proceeds  as  follows: 


T  a  T  A  T  A 

x.  Q  x.  ♦  2  u.  M  x.  ♦  u.  R  u. 

—l  —i  —l  —l  —l  —l 


f  i+1 

=  I  l(*£i  *  IV)T  Qi^Xj  +  IV)  +  2  u^ 

*i 

T 

♦  Uj  R  Uj]  dt 

f^i+l  T  t  f^i+1  T  T 

*  (£i  *  Q  *  *i)dt  ♦  2  (Hi  r  Q  *  £i 

*i  n 

t  f*i+l  T  T  T 

+  u.M  4>  x.  )dt  ♦  [u.  R  u.  ♦  u.  r*Qr  u. 

—1  — 1  J  — 1  —1  —x  —1 

t. 

l 

♦  2  (u.TMr  u.}]  dt 


Because  x^  and  u^  are  not  functions  of  the  integration  variable,  t, 
they  may  be  taken  outside  the  integrals.  Matching  like  variables 
from  both  sides  of  the  equation  we  have  the  desired  relation: 
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APPENDIX  B 


DESCRIPTION  OF  APL  FUNCTIONS  FOR 
GENERATING  OPTIMAL  GAINS  AND  TIME  HISTORIES 


The  equations  from  modern  control  theory  (given  in  Chapter  2) 
are  coded  in  APL  functions.  The  functions  are  listed  in  Table  B-l. 

A  short  description  of  each  function  is  given  below: 


FCLOOP  -  Calculates  the  closed-loop  F  matrix  for  a 

second-order  system  given  4>,  T,  and  the  gains, 
k.  F(;l  is  returned  in  variable  FCL.  Uses 
Eq.  46. 

FCL4  -  Calculates  the  closed- loop  F  matrix  for  a 
fourth-order  system  given  the  fourth-order 
$  and  T,  and  the  second-order  gains,  k. 

Fcl  is  returned  in  variable  FCL. 

FOLLOW  -  Calculates  continuous-time  Q,  M,  and  R  for 
implicit  model-following  given  the  matrix 
F  of  the  model,  matrices  F  and  G  of  the  VRA 
and  a  weighting  matrix,  W.  Uses  Eqs.  57-59. 

GAIN  -  Calculates  optimal  gains  using  Eq.  30. 

Gains  are  passed  to  variable  k.  The  di¬ 
mension  of  k  is  changed  in  k4  for  use  in 
FCL4. 

GAMMA  -  Calculates  the  discrete  control  effects 

matrix,  T,  for  a  second-order  system.  It 
is  a  monadic  operator  that  uses  the  second- 
order  4>(STM).  Variable  GAM  holds  the  re¬ 
sult.  Uses  Eq.  12. 

GAMMA4  -  Calculates  fourth-order  T.  It  is  a  monadic 
operator  that  uses  the  fourth-order  4>(STM4). 

GENSTM  -  Generates  100  $  matrices  using  time  inter¬ 
vals  in  1  percent  increments  of  the  sam¬ 
pling  interval.  These  matrices  are  for 
use  in  function  QMR  for  calculating  dis¬ 
crete  weighting  matrices. 


r 


IDENT  -  Forms  an  identity  matrix.  A  monadic 

operator  that  uses  the  matrix  dimension 
as  the  argument.  Result  returned  in 
variable  I. 

QMR  -  Qalculates^discrete  weighting  matrices 
Q,  M,  and  R  from  continuous  Q  and  R. 
Uses  simple  Euler  integration  in  100 
steps.  Uses  Eqs.  26-28. 

QMRF  -  Calculates  discrete  weighting  matrices 

Q,  and  6  from  continuous  Q,  M,  and 

R.  Uses  simple  Euler  integration  in 
100  steps.  Uses  Eqs.  60-62. 

RICCATI-  Solves  Riccati  equation  (Eq.  29). 

Result  returned  in  variable  P. 

SONETWO-  Calculates  relation  between  Ay^  and 
Ax*:Si2.  It  uses  the  result  from 
STWOTWO.  Uses  Eq.  18. 

STMDYN  -  Calculates  $  given  sampling  interval 
and  F  matrix.  Result  is  returned  to 
variable  specified  by  using  the  "back 
arrow."  Uses  Eq,  11. 

STWOTWO-  Calculates  relation  between  Ay^  and  Au*: 
S22-  Uses  Eq.  19,  ~t  ~ 


f  V 


I 


t  SYSMODEL- Generates  a  time  history  for  the  second- 

order  system.  The  function  will  ask  if 
a  printed  table  is  desired  (answer  yes 
or  no).  The  time  history  will  run  for 
5  sec.  unless  'BREAK'  is  depressed.  The 
time  history  is  stored  in  vectors  AQ, 
t  AALPHA,  TIME,  and  DELU.  Vector  AQ  con¬ 

tains  elements  Axj^.  AALPHA  contains 
Ax2k,  TIME  contains  the  time  instances 
spaced  by  the  sampling  interval,  and 
DELU  contains  Auj^.  These  vectors  are 
used  in  the  10  LINPLOT  functions  to  plot 
i  the  time  histories.  Uses  Eqs.  10  and  22. 


SYSMOD4-Same  as  SYSMODEL  except  4th-order  time  history 
is  generated  using  the  2nd-order  gains.  Addi¬ 
tional  vectors  for  plotting  are  Av  and  AGAM. 

Vector  Av  contains  Ax»  and  AGAM  contains  Axa,  . 
»  -^k  k 
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There  are  two  separate  methods  for  finding  the  optimal 
gains  using  these  functions.  The  first  method  involves  changing 
the  weighting  matrices  in  the  cost  functional  to  obtain  the  desired 
response.  The  sequence  for  using  the  APL  functions  in  this  manner 
is  given  in  Fig.  B-l.  For  example,  say  the  gains  are  to  be  found 
for  a  pitch  rate  controller  using  the  second-order  model  given  in 
Eq.  39.  It  can  be  specified  as  a  pitch  rate  controller  by  letting 


Hx  =  [1  0]  and  Hu  =  0 


The  weightings  for  this  example  will  be 


25  0 

Q  *  ,  fi  =  33 

0  25 


and  a  sampling  interval  of  0.1  sec  is  used. 


The  public  workspace  3  EIGENVAL  is  copied  for  use  in 
calculating  eigenvalues  from  F^-  Following  the  sequence  of 
Fig.  B-l: 
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Figure  B-l.  Sequence  for  Calculating  Gains  by  Changing 
Weighting  Matrices. 


Figure  B*l. 


continued 


The  second  method  for  calculating  gains  uses  implicit 
model-following.  For  instance,  say  the  VRA  is  to  match  the  short- 
period  modal  characteristics  of  an  aircraft  with  the  reduced  second- 


order  model 


«■[•■:  xi  •  n  - 


A  weighting  matrix,  W,  should  initially  be  set  to  the  identity  matrix. 

s 

Following  the  sequence  given  in  Fig.  B-2: 
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Figure  B-2.  Sequence  for  Calculating  Gains  Using 
Implicit  Model*Following. 
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RESEARCH  SYSTEMS 

C.l  VARIABLE-RESPONSE  RESEARCH  AIRCRAFT  (VRA) 

The  VRA  is  a  highly  modified  Navion  equipped  with  inertial, 
air  data,  and  navigation  sensors,  as  well  as  six  independent  force 
and  moment  controls.  The  VRA,  shown  in  Fig.  C-l,  has  been  used  to 
conduct  a  broad  range  of  experiments  in  aircraft  flying  qualities, 
human  factors  and  control  in  the  past.  Tht  aircraft  has  played  a 
major  role  in  establishing  current  military  and  civil  flying  qual¬ 
ities  criteria,  and  with  the  addition  of  the  Micro-DFCS,  the  VRA  is 
equipped  to  expand  this  type  of  research,  as  well  as  to  investigate 
advanced  digital  control  concepts  (Ref.  23) . 

Independent  control  of  three  forces  and  three  moments  is 
provided  by  commands  to  the  elevator,  ailerons,  rudder,  throttle, 
direct-lift  flaps,  and  side-force  panels  (Fig.  C-2).  The  control 
surfaces  are  driven  by  hydraulic  servos  originally  fitted  to  the 
B-58  aircraft.  The  modified  VRA  units  incorporate  solenoid-actuated 
valves  with  force-override  features  for  quick  disengagement. 
Characteristics  of  the  control  effectors  are  summarized  in  Table 
C-l.  Surface  rate  limits  are  seen  to  range  from  60  to  110  deg/sec. 
Bandwidths  are  given  for  flat  response  and  6  db  attenuation  (in 
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Figure  C-2.  Major  Components  of  VRA. 


parentheses),  except  that  thrust  bandwidth  is  specified  by  the  fre¬ 
quency  for  3  db  attenuation.  The  aircraft's  normal  operating  speed 
range  is  65  to  120  kt;  maximum  specific  forces  and  moments  ("control 
power")  are  given  for  70  kt  airspeed.  At  IAS  =  105  kt,  maximum 
direct  lift  and  side-force  accelerations  are  lg  and  0.5g,  respectively. 

The  sensors  used  for  most  flight  testing  include  angular 
rate  gyros  and  linear  accelerometers  for  all  three  axes,  vertical  and 
heading  gyros,  dual  angle-of-attack  and  sideslip-angle  vanes,  radar 
altimeter,  indicated  airspeed,  control  surfacs  positions,  and  cockpit 
control  positions.  Several  other  signals  (e.  j.,  air  temperature, 
barometric  altimeter,  altitude  rate,  and  TALAR  microwave  landing 
system  signals)  are  available  for  system  feedback  or  telemetry  re¬ 
cording.  The  present  telemetry  system  allows  42  data  channels 
to  be  multiplexed  and  transmitted  to  the  FRL  ground  station  described 
below. 


TABLE  C-l 

VRA  Control  Characteristics 


* 

Control 

Displacement 
Limit,  deg 

Rate  Limit, 
deg/sec 

Bandwidth, 

Hz 

Maximum  Specific 
Force  or  Moment 
(IAS  =  70kt) 

Roll 

30. 

70. 

5  (10) 

4.1  rad/sec^ 

» 

Pitch 

-30. 

♦  15 

70. 

5  (10) 

2 

4.4  rad/sec 

Yaw 

15. 

70. 

5  (10) 

1.9  rad/sec^ 

Thrust 

- 

- 

0.6 

0.1  g 

1 

Side 

Force 

35. 

60. 

2  (3) 

0.25  g 

Normal 

Force 

30. 

110. 

2  (3) 

0.5  g 

I 


1 


t” 

l 


f 


... 


)**'•  - 
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The  aircraft  is  flown  by  a  two-man  crew  during  all  research. 
This  provides  a  number  of  advantages  in  comparison  to  single-pilot 
operation  from  the  standpoint  of  flight  safety  and  experimental 
efficiency.  The  instrument  panel  and  controls  are  shown  in  Fig.  C-3. 
The  conventional  mechanical  aircraft  system  is  flown  by  the  safety 
pilot  in  the  left  seat,  while  the  fly-by-wire  aircraft  system  used 
for  research  is  flown  by  the  evaluation  pilot  seated  at  the  right. 
This  system  includes  the  Micro-DFCS  and  redundant  aileron,  elevator, 
and  side-force  actuators  for  protection  against  system  failures. 

The  evaluation  pilot's  station  is  tailored  to  the  experiment;  for  the 
Micro-DFCS  program,  this  station  includes  a  center  control  stick, 
thumb  switches  for  trim  and  direct  force  modes,  rudder  pedals, 


Figure  C-3.  Instrument  Panel  and  Controls  of  VRA. 
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sideslip  and  side-force-panel  meters,  and  conventional  instruments. 

The  safety  pilot  is  the  in-flight  test  conductor,  monitoring 
systems  and  adjusting  all  experimental  parameters.  He  has  several 
electrical  and  hydraulic  mechanisms  for  disengaging  the  Micro-DFCS 
and  the  variable-response  system  in  the  event  of  a  malfunction,  as 
well  as  an  "automatic  go-around"  abort  mode  which  makes  safe  experi¬ 
mentation  through  touchdown  possible.  The  abort  mode  commands  a 
20-deg  flap  setting  and  climb  power  when  activated;  at  70-kt  (36  m/s) 
airspeed  on  a  6-deg  glideslope,  an  up- flap  "hardover"  failure  can  be 
corrected  and  climbout  can  be  initiated  with  a  maximum  altitude  loss 
of  10  ft  (3  m) . 

C. 2  EXPERIMENTAL  FACILITIES 

The  VRA  is  operated  from  the  flight  test  facility  at 
Princeton  University's  James  Forrestal  Campus.  The  facility  includes 
the  FRL  hangar,  laboratories,  and  shops,  plus  a  3000-ft  Basic  Utility 
II  runway.  TALAR  3  and  4  fixed-beam  microwave  landing  systems  (MLS) 
furnish  precision  approach-path  guidance. 

The  ground  station  (Fig.  C-4)  at  the  FRL  is  used  to  receive, 
record,  and  analyze  the  telemetered  data  from  the  VRA.  It  includes 
a  Honeywell  seven  channel  tape  recorder,  an  FM  or  AM  receiver  pre¬ 
sently  operating  at  1458  MHz  in  the  FM  mode,  a  telemetry  demulti¬ 
plexer  with  five  translators,  an  EAI  TR-48  analog  computer,  a  radio 
telephone,  and  a  six  channel  paper  strip  chart  recorder.  The  PDM 
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telemetry  system  provides  42  data  channels,  each  sampled  at  a  rate 
of  20  sps.  The  telemetry  data  from  the  receiver  can  be  recorded 
on  tape  and  demultiplexed  five  channels  at  a  time  for  plotting  on 
the  strip  chart  recorder.  The  analog  computer  scales  and  buffers 
all  input  channels  from  the  translators  to  the  strip  chart  recorder. 


Figure  C-4.  FRL  Ground  Station. 


In  addition  to  the  analog  computer's  function  in  analyzing 
telemetered  data,  it  is  also  used  for  ground-based  simulations  of 
the  VRA  and  other  dynamic  systems.  In  the  Micro-DFCS  study,,  it  is 
used  to  model  the  VRA  for  ground  testing  the  microcomputer  software. 
The  schematic  and  corresponding  potentiometer  settings  are  presented 
in  Fig.  C-5  and  Table  C-2,  respectively. 
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Figure  C-5.  Schematic  of  Simulated  VRA  at  105  Kts 

With  Scaling  for  Microcomputer  Interface. 


TABLE  C-2 


COMPUTER  POTENTIOMETER  SETTINGS 
FOR  VRA  LONGITUDINAL  MODEL  (105  KIAS) 


Pot 

Parameter 

Scaling 

Setting 

00 

m<Se 

M6e/100 

.125 

01 

Ma 

V10 

.835 

02 

Me 

Mg/10 

.208 

03 

Md 

.910 

05 

%/v 

.125 

06 

La/V 

La/v/iO 

.200 

07 

LV/V 

100  Vv 

.200 

08 

Step  Gust 

Scaling 

.100 

10 

D  -T 

V  V 

.073 

11 

Da-g 

Da-g/l00 

.060 

12 

g 

g/100 

.322 

13 

vscale 

100 

57.3x5 

.349 

25 

V/g 

v/  *  2x10 

V/8  *573 

.190 

30 

scaling  error 

“MSE/io 

.160 

31 

1 

ACONV 

1 

.345 

32 

1 

.144 

ucOnv 

UCONVxlO 

33 

1 

QCONV 

.225 

47 

1 

nconv 

.188 

55 

Stick  Range 

.625 

56 

Stick  Scaling 

.592 

Nominal  Scaling:  1  volt  *  1°,  l°/sec,  5  ft/sec,  .05  G's 
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THE  MICRO-DFCS  HARDWARE 

D.l  DESCRIPTION  OF  MICROCOMPUTER  COMPONENTS 

In  addition  to  the  brief  description  of  the  microcomputer 
equipment  given  in  Chapter  3,  the  following  information  is  provided. 
The  CPU  board  has  22  parallel  lines  of  input  and  output  (I/O),  4.5 K 
of  memory  and  one  hardware  interval  timer  which  may  be  wired  to  an 
irterrupt  line  of  the  8085.  The  timer's  time  interval  may  be  changed 
under  software  control.  The  mathematics  unit  does  16-  and  32-bit 
fixed-point  and  32-bit  floating-point  operations,  where  the  typical 
time  needed  to  do  one  floating-point  operation  (including  the  time 
to  pass  arguments  to  the  mathematics  unit)  is  137  usee.  The  mathe¬ 
matics  unit  provides  the  computational  speed  necessary  to  calculate 
the  control  law  within  the  sampling  interval.  The  main  memory  con¬ 
tains  16K  words  of  RAM  and  also  provides  an  interface  for  the  CDU. 

The  battery  back-up  memory  allows  programs  and  data  to  be  stored 
with  the  power  off  for  up  to  96  hours.  This  board  is  used  to  trans¬ 
fer  the  microcomputer  programs  from  the  software  development  system 
to  the  FCCU  mounted  in  the  VRA.  The  analog  I/O  boards  provide  16 
(expandable  to  32)  input  channels  and  6  output  channels  each  with 
12-bits  of  resolution.  The  A/D  has  a  conversion  rate  of  28KHz  and 
can  be  interrupt-  or  software-driven.  The  CDU  provides  input  for 
all  ASCII  characters  and  can  display  two  lines  of  12  characters  each. 


It  is  functionally  equivalent  to  a  conventional  1200-baud  keyboard/ 
display  terminal,  although  its  display  is  limited,  and  multiple  key 
strokes  are  required  to  enter  most  characters. 

The  FCCU  consists  of  the  six  computer  circuit  boards  and 
two  4-board  cardcages  into  which  the  boards  plug.  It  is  housed  in 
an  RF-shielded,  shock-mounted  aluminum  box.  The  mount  is  located 
in  the  cockpit  behind  the  two  pilots.  The  FCCU  therefore  experiences 
the  same  environment  the  pilots  do,  and  no  special  provisions  are 
made  in  using  these  commercially  available  microcomputer  components. 
The  Micro-DFCS  components  installed  in  the  VRA,  including  the  power 
supplies,  cost  approximately  $6900  when  purchased  in  February,  1978. 
There  are  four  cables  that  interface  the  FCCU  with  the  rest  of  the 
VRA's  FBW  system.  One  cable  supplies  the  +5  and  +12  volts  necessary 
to  operate  the  Intel  SBC  boards.  These  voltages  are  obtained  by 
regulating  the  VRA's  primary  28vdc.  A  second  connector  links  the 
CDU  to  the  FCCU.  A  two-position  switch  and  program  monitor  light 
mounted  on  the  pilot's  instrument  panel  each  have  connections  made 
through  a  third  cable.  One  position  of  the  switch  resets  the  com¬ 
puter's  program  counter,  and  the  other  position  gives  the  CPU  program 
control.  The  light  is  part  of  an  error  detection  system.  When  the 
light  flashes  at  a  steady  rate,  it  indicates  that  the  control  law 
is  being  calculated  every  sampling  interval.  The  fourth  cable  sup¬ 
plies  the  inputs  from  the  motion  sensors  and  pilot  controls,  and  the 
outputs  to  the  control  surface  actuators.  The  voltage  range  of  these 
signals  (+10v)  is  compatible  with  that  used  by  the  A/D  and  D/A  con¬ 
verters  of  the  microcomputer. 


The  status  of  the  Micro-DFCS  is  output  on  one  of  the  D/A 
channels  and  this  information  is  sent  on  a  telemetry  channel.  This 
status  channel  indicates  what  control  mode  is  running  (if  any)  and 
also  flags  any  errors  found  by  the  two  software  error  detection 
methods  described  in  Section  3.2.  The  telemetered  data  can  be  re¬ 
corded  and  plotted  on  a  strip  chart  recorder. 

The  components  of  the  Micro-DFCS  and  development  system 
are  pictured  in  Fig.  D-l.  The  FCCU  and  CDU  are  at  the  upper  left, 
resting  on  the  ground  chassis  and  power  supply.  The  keyboard-CRT 
terminal,  acoustic  coupler,  and  telephone  extension  are  to  the 
right  of  the  chassis.  A  keyboard-printer  unit  (not  shown  in  Fig.  D-l) 
is  used  to  obtain  copies  of  program  listings  and  other  computer  output. 


Figure  D-l.  Components  of  the  Micro-DFCS 

and  Software  Development  System. 
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D.2  DIRECTIONS  FOR  USING  CONTROL  DISPLAY  UNIT 


The  Control  Display  Unit  (CDU) ,  shown  in  Fig.  D-2,  allows 
the  pilot  to  monitor  the  Micro-DFCS  and  input  commands.  The  dis¬ 
play  is  two  lines  of  12  characters  each.  The  keyboard  allows  all 
128  ASCII  characters  to  be  entered;  most  by  multiple  key  strokes. 


FUNCTION 

Miniature  hand  held  computer  terminal  enables  bi-directional  communication,  using 
ASCII  codes  in  a  bit  senal,  asynchronous  format. 


f  DISPLAY 

Computer  output  (or  input/output)  is  displayed 
here  in  one  of  several  formats 


□adJdl' 


SHIFT  KEYS 

Operated  by  finger*  of 
the  left  hand  to  change 
interpretation  of 
KEYPAD  keys 


6—  ft  SIGNAL 
and  POWER  Coed 


KEYPAD 

Twenty  push  keys  opeiated  by 
fingers  of  the  tight  hand  Each 
key  has  several  character 
interpretations  as  conttolled  by 
the  three  SHIFT  KEYS  and  a 
CASE/CONTROL  key  Enables 
generation  of  all  1 28  ASCII 
characters  plus  “break" 


Figure  D-2.  Control  Display  Unit  (CDU) 


There  are  three  levels 
System  Monitor  level  is  entered 
on  the  instrument  panel  (marked 


of  entry  on  the  CDU  (Fig.  D-3).  The 
by  flipping  the  two-position  switch 
RESET/RUN)  to  the  RUN  position.  The 
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Figure  D-3.  Levels  of  Entry  on  CDU. 

CDU  display  will  show  "80/05  MON"  indicating  the  first  entry  level. 
Any  of  the  normal  80/05  System  Monitor  commands  outlined  in  the 
Intel  SBC  80P05  User's  Guide  can  be  used  at  this  level. 

The  second  level  of  entry  may  be  entered  by  jumping  to 
location  8000H.  This  is  accomplished  by  pressing  the  keys: 

G,  8,  0,  0,  0,  CR 

where  CR  is  the  CDU's  "carriage  return".  This  starts  execution  of 
CAS-1  and  enters  the  CAS-1  Command  Recognition  level.  Twenty-six 
different  commands  may  be  entered  at  this  level,  each  selected  by 
pressing  a  key  with  a  letter  of  the  alphabet.  Only  four  commands 
are  presently  defined  in  CAS-1:  HALT,  INITIALIZE,  MODE  CHANGE,  and 
STEP  INPUT. 
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The  third  level  sets  parameters  necessary  to  complete  the 
command  entered  at  the  second  level.  For  example,  if  MODE  CHANGE 
was  entered  by  pressing  "M",  then  the  third  entry  level  is  auto¬ 
matically  set  and  the  next  two  key  strokes  will  specify  which  mode 
is  selected  (a  number  00-19).  Once  the  mode  number  is  entered,  then 
the  CDU  again  operates  at  the  second  entry  level  waiting  for  another 
command.  STEP  INPUT  also  requires  the  third  entry  level  after  the 
"S"  is  pressed.  The  analog  channel  to  be  stepped  (00-15)  is  entered 
followed  by  the  voltage  this  channel  is  to  be  stepped  (in  the  format 
+DD.DD,  where  D  represents  some  base  ten  number). 

The  upper  line  of  the  display  shows  the  mode  CAS-1  is  oper¬ 
ating  in  at  all  times.  In  the  case  where  CAS-1  is  first  entered, 
the  mode  is  Initialization  and  the  top  display  line  reads  "8085 
READY".  If  the  mode  is  changed  by  using  the  MODE  CHANGE  command, 
the  new  mode  is  displayed  on  the  upper  line  (e.g.  "OPEN  LOOP", 

PITCH  RATE",  or  "NORMAL  ACC").  The  entries  made  at  the  second  level 
are  displayed  on  the  bottom  line.  If  an  entry  error  is  made  at  any 
time  while  using  the  CDU,  "ENTRY  ERROR"  will  be  displayed  on  the 
bottom  line  for  ten  sampling  intervals. 

To  reenter  the  System  Monitor  level,  the  switch  on  the 
instrument  panel  is  placed  in  the  RESET  position  and  then  back  to 
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THE  MICRO-DFCS  SOFTWARE 

E.l  DESCRIPTION  OF  CAS-1  ROUTINES 

The  routines  of  CAS-1  are  divided  into  three  categories, 
depending  on  their  function:  Executive,  Utility  and  Flight  Control. 
Each  of  these  routines  has  detailed  documentation  included  in  the 
source  listing  consisting  of  line  comments  and  a  header.  The  header 
explains  the  purpose  of  the  routine,  how  variables  should  be  passed 
to  it,  and  what  parameters  of  the  microcomputer  are  effected  by  the 
routine's  execution. 

A  brief  explanation  of  the  major  elements  of  CAS-1  is  given 

below. 

EXECUTIVE  ROUTINES  -  Initialization,  CDU  interface,  and 
one  of  the  three  error  detection  methods  are  contained  in  this  sec¬ 
tion  (Fig.  E-l).  INITIALIZATION  involves  setting  up  the  mathematics 
unit,  the  analog  board,  the  parallel  I/O  ports,  the  serial  port, 
the  hardware  interrupt  timer,  the  registers  to  be  used  as  counters 
and  flags,  and  clearing  storage  area  in  RAM.  CDU  INTERFACE  routine 
allows  the  user  to  set  the  operation  of  CAS-1  by  inputs  from  the 
keyboard.  MEMORY  CHECK  is  one  of  two  software  error  detection 
methods  employed  in  the  program.  Every  50  samples,  the  entire  con¬ 
tents  of  CAS-1  is  summed  in  an  8-bit  register  (ignoring  overflow) 
and  compared  to  the  known  sum.  If  the  two  differ,  some  part  of 
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the  program  has  changed;  this  condition  is  displayed  on  the  COU, 
and  sent  via  telemetry  to  the  ground  station. 

UTILITY  ROUTINES  -  The  routines  CLEAR  LINE,  CONSOLE  OUTPUT, 
NUMERIC  INPUT,  ENTRY  ERROR,  and  SERIAL  OUTPUT  are  used  for  displaying 
and  entering  data  on  the  CDU.  BLINK  is  part  of  the  hardware  error 
detection  method  that  flashes  a  light  on  the  safety  pilot's  instru¬ 
ment  panel  to  indicate  that  the  flight  control  function  of  CAS-1  is 
running.  BLINK  is  called  once  in  every  flight  control  interrupt 
service  routine.  ANALOG  TO  DIGITAL  CONVERSION  selects  the  proper 
analog  input  channel,  initiates  the  conversion,  and  stores  the  re¬ 
sults.  COUNT-UP  DISPLAY  generates  an  increasing  sequence  of  numbers 
1-9  (reset  each  10  sec)  on  the  bottom  line  of  the  CDU.  This  is  to 
indicate  that  the  hardware  interrupt  (used  for  timing  purposes)  is 
working,  the  program  has  initialized  correctly,  and  the  D/A  channels 
are  operational.  An  increasing  sequence  of  voltage  outputs  from  a 
D/A  channel  is  telemetered  to  the  ground  station  to  indicate  in 
what  mode  the  Micro-DFCS  is  operating.  INTERRUPT  COUNT  increments 
a  counter  every  time  an  interrupt  service  routine  is  executed.  This 
counter  is  used  for  timing  purposes  in  routines  such  as  MEMORY  CHECK 
(which  must  execute  every  50  samples) .  LIMIT  ANALOG  OUTPUT  prevents 
the  analog  output  channels  from  instantly  switching  from  +10  volts 
to  -10  (or  vice  versa).  This  condition  occurs  when  the  control  law 
calculates  a  value  for  control  surface  deflection  which  corresponds 
to  a  voltage  greater  than  10  volts  or  less  than  -10  volts.  It  is 
because  of  the  method  used  for  converting  floating  point  numbers  to 
fixed  point  format.  This  routine  was  not  added  until  after  flight 


testing  had  begun.  MODE  CHANGE  allows  the  user  to  select  1  of  20 
possible  flight  control  inodes.  Presently  only  six  modes  are  de¬ 
fined;  four  pitch  rate  command  modes  with  different  feedback  gains, 
one  normal  acceleration  command  mode,  and  one  direct  mode  where 
longitudinal  stick  is  fed  directly  to  the  elevator  through  a  stick 
gearing  factor.  MATH  UNIT  DRIVER  loads  the  high-speed  mathematics 
unit  with  the  two  arguments  to  be  operated  on,  initiates  the  mathe¬ 
matical  operation,  and  stores  the  results  in  RAM.  It  is  very  gen¬ 
eral  and  slow  (normally  230  usee).  Therefore,  it  should  not  be 
used  for  time  critical  functions;  i.e.  not  in  flight  control  inter¬ 
rupt  service  routines.  MATH  ERROR  PROCESSOR  is  the  second  software 
error  detection  method.  It  is  used  to  check  for  errors  that  occur 
in  the  mathematics  unit  (such  as  divide  by  zero  or  an  overflow  con¬ 
dition)  .  The  type  of  error  that  has  occurred  and  the  number  of 
mathematics  unit  uses  since  the  last  interrupt  service  routine  be¬ 
gan  execution  is  displayed  on  the  CDU.  The  type  of  error  is  also 
sent  to  the  ground  station  via  telemetry.  This  information  should 
allow  the  user  to  find  which  mathematical  operation  is  causing  the 
problem,  should  one  exist.  CALIBRATED  STEP  INPUT  allows  the  pilot 
to  apply  a  step  input  on  any  one  of  the  analog  input  channels. 

After  entering  the  desired  step  value,  the  input  is  initiated  by 
depressing  the  CDU's  "carriage  return"  key.  The  step  input  is 
nulled  by  depressing  any  CDU  key. 

FLIGHT  CONTROL  ROUTINES  -  These  are  of  two  types:  control 
set-up  routines  and  the  timed  interrupt  service  routines.  The  con¬ 
trol  set-up  routines  are  called  from  MODE  CHANGE  and  set  the 


parameters  necessary  for  the  interrupt  service  routines  to  operate. 

A  set-up  routine  will  store  the  starting  location  of  a  service  rou¬ 
tine  in  the  timed  interrupt  branch  point,  and  it  will  set  the  op¬ 
timal  gains  for  the  control  law.  A  set-up  routine  will  also  store 
the  nominal  values  of  all  analog  inputs  for  use  in  calculating  per¬ 
turbation  values  by  the  service  routine.  The  microcomputer  status 
channel  is  set  to  indicate  what  flight  control  mode  the  computer  is 
executing. 

The  flight  control  interrupt  service  routines  are  executed 
on  every  timed  interrupt.  They  contain  the  control  logic  that  cal¬ 
culates  the  control  inputs,  Au,  for  every  sampling  interval.  Once 
a  set-up  routine  has  set  the  timed  interrupt  branch  point  to  exe¬ 
cute  a  particular  service  routine  on  every  timed  interrupt,  then 
that  service  routine  has  the  highest  priority  over  all  other  routines 
and  will  be  executed  on  every  timed  interrupt,  regardless  of  where 
the  program  is  running  at  the  time. 


The  DIRECT  INTERRUPT  SERVICE  ROUTINE  takes  in  commands 
from  longitudinal  stick,  adds  the  step  bias,  multiplies  this  input 
by  a  stick  gearing  factor,  and  sends  this  value  to  the  elevator 
actuator.  The  stick  gearing  used  for  flight  tests  was  three  degrees 
elevator  deflection  per  inch  of  longitudinal  stick,  6g,  deflection. 
At  105  KIAS,  this  corresponds  to  a  steady-state  Aq  of  5.3  deg/sec 
per  inch  or  a  steady-state  An2  of  0.52g  per  inch  of  6S. 

The  PITCH  RATE  INTERRUPT  SERVICE  ROUTINE  calculates  an  ele¬ 
vator  command,  6g,  for  measured  values  of  6S,  o,  and  q.  The  inputs 
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must  be  formatted  and  scaled  before  the  control  law  can  begin  its 
calculation  (Fig.  E-2).  First  the  analog  voltage  signal  for  a 
particular  input  is  converted  to  a  12-bit  fixed-point  binary  num¬ 
ber,  which  is  then  converted  to  a  32 -bit  floating-point  represen¬ 
tation.  The  nominal  value  for  the  channel  (stored  by  the  control 
set-up  routine)  is  subtracted  from  the  current  input  value  to  form 
the  perturbation  value.  At  this  point,  if  a  step  bias  value  has 
been  specified  for  a  particular  input  channel  by  CALIBRATED  STEP 
INPUT,  that  value  is  added  to  the  input.  The  control  law  must  be 
calculated  with  the  inputs  scaled  to  the  proper  units.  A  scaling 
factor  that  converts  units  of  voltage  to  the  proper  units  of  motion 
is  the  last  step  before  actual  calculation  of  Au.  In  the  case  of 
the  pilot's  inputs,  this  scaling  factor  can  be  thought  of  as  a  stick 
gearing.  For  Aa,  there  is  one  more  scaling  to  be  done  that  takes 
into  account  the  difference  between  the  a  vane  deflection  and  the 
aircraft's  actual  a.  The  output  from  the  control  law  is  converted 
from  units  of  radians  to  volts,  converted  to  fixed-point  notation, 
then  converted  to  an  analog  voltage,  and  sent  to  6g. 

The  control  law  is  calculated  using  Eq.  22  where  the  mathe¬ 
matical  operations  are  performed  by  calling  MATH  UNIT  DRIVER.  This 
mathematics  routine  is  actually  too  slow  and  general  for  calculating 
the  control  law.  The  PITCH  RATE  INTERRUPT  SERVICE  ROUTINE  takes, 
on  the  average,  19  msec  to  execute  including  all  formating  and 
scaling.  The  normal  acceleration  command  version  of  this  routine 
is  written  not  using  MATH  UNIT  DRIVER  in  order  to  decrease  execu¬ 
tion  time. 
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The  NORMAL  ACCELERATION  INTERRUPT  SERVICE  ROUTINE  calcu¬ 
lates  a  $£  command  for  measured  values  of  6S>  nz,  and  q.  The  task 
of  this  routine  is  similar  to  the  pitch  rate  service  routine,  but 
the  method  of  making  calculations  is  fundamentally  different  in 
two  ways.  First,  the  MATH  UNIT  DRIVER  is  not  used  to  perform  the 
calculations.  Instead, this  service  routine  sets  up  the  mathematics 
unit  and  initiates  the  mathematical  operations  itself.  While  the 
mathematics  unit  is  doing  the  calculation,  the  service  routine  gets 
the  next  calculation  ready.  Also,  only  one  argument  is  loaded  to 
the  mathenatics  unit  per  calculation  because  the  answer  from  the 
last  calculation  is  always  used  as  the  second  argument.  By  compar¬ 
ison,  MATH  UNIT  DRIVER  always  loads  two  arguments  and  stores  the 
answer  regardless  of  the  next  operation  to  be  performed.  By  not 
using  MATH  UNIT  DRIVER,  this  service  routine  takes  more  storage 
but  each  calculation  is  speeded  up  by  a  factor  1.7  (137  ysec  per 
calculation).  The  second  difference  in  calculations  comes  from 
using  a  reduced  form  of  Eq.  22  for  the  control  law.  The  equation 
can  be  rearranged  to  the  form 

=  (S22  +  CSl2^  Zdk  ‘  C  AXfc  (74) 

For  the  two-state  feedback  and  single  pilot  input  CAS  studied  here, 
Eq.  74  breaks  down  further  to 

A6Ek  “  cj  A6S  -  c2  Aq  -  Cj  A  nz  (75) 
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which  only  takes  five  mathematical  operations  to  calculate  compared 
to  nine  for  Eq.  22.  By  not  using  MATH  UNIT  DRIVER  and  using  Eq.  75 
for  the  control  law,  execution  time  of  the  NORMAL  ACCELERATION 
INTERRUPT  SERVICE  ROUTINE  is  reduced  to  6  msec.  PITCH  RATE  INTERRUPT 
SERVICE  ROUTINE  could  also  be  rewritten  to  run  at  this  speed  because 
it  has  the  same  number  of  calculations  as  the  normal  acceleration 
routine.  The  scaling  factors  for  the  inputs  and  outputs  can  be  in¬ 
corporated  in  the  gains  of  Eq.  75  to  speed  up  the  service  routines 
even  more.  This  was  not  done  in  CAS-1  to  help  isolate  any  problems 
that  might  have  occurred  in  its  initial  use. 

E.2  SOFTWARE  DEVELOPMENT 

Each  routine  fox  the  Micro-DFCS  software  was  developed 
using  the  method  shown  in  Fig.  E-3.  The  routines  are  written  in  the 
assembly  language  for  the  Intel  8085  CPU.  This  language  consists 
of  80  different  mnemonics  that  represent  the  machine  operation  codes. 
The  assembly  language  is  cross-assembled  to  the  machine  language  and 
loaded  to  the  microcomputer  for  debugging.  Initially,  the  routines 
are  tested  separately,  i.e.,  without  the  other  program  modules. 

This  procedure  helps  identify  problems  that  may  occur  when  the  rou¬ 
tines  are  combined  to  operate  in  the  program. 

E.3  ADDING  FLIGHT  CONTROL  ROUTINES  TO  CAS-1 

New  or  modified  Flight  Control  routines  can  easily  be  added 
to  CAS-1.  The  procedure  consists  of  l)  defining  a  new  mode  in  the 
Mode  Branch  Table  of  CAS-1,  2)  writing  an  appropriate  mode  set-up 
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routine,  and  3)  writing  the  flight  control  interrupt  service  rou¬ 
tine  in  a  form  that  makes  it  compatible  with  the  operation  of  CAS-1. 
Step  one  is  the  only  modification  necessary  to  the  coding  of  CAS-1. 
It  consists  of  replacing  the  operand  "BADEN"  at  the  correct  lo¬ 
cation  in  the  Mode  Branch  Table  with  the  label  for  the  new  mode 
set-up  routine. 

In  the  second  step,  the  mode  set-up  routine  is  written  to 
set  the  timed  interrupt  branch  point  and  store  parameters  for  the 
interrupt  service  routine.  The  branch  point  for  the  timed  interrupt 
starts  at  3FECH  in  RAM.  A  typical  mode  set-up  routine  is  structured 
as  in  Fig.  E-4.  All  registers  are  saved  to  insure  the  routine  does 
not  change  any  parameters  for  CAS-1.  The  interrupt  is  disabled  so 
the  interrupt  service  routine  does  not  execute  before  the  gains  and 
other  parameters  can  be  changed.  The  mode  type  (e.g.  OPEN  LOOP, 
PITCH  RATE,  etc.)  is  displayed  on  the  bottom  line  of  the  display  and 
a  corresponding  voltage  level  is  set  on  the  analog  status  channel. 
After  storage  is  set  for  the  service  routine,  all  registers  are  re¬ 
stored  and  the  interrupt  is  enabled  to  allow  the  newly  specified 
service  routine  to  execute  on  the  occurence  of  each  timed  interrupt. 

The  third  step  of  the  procedure  is  to  structure  the  flight 
control  interrupt  service  routine  as  shown  in  Fig.  E-5.  The  inter¬ 
rupt  is  disabled  to  guarantee  the  routine  completes  before  the  next 
service  routine  starts.  All  registers  are  saved  because  the  timed 
interrupt  may  break  the  program  flow  of  CAS-1  at  any  point,  and 
therefore  the  status  of  the  program  at  the  break  must  be  saved. 
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Output  port  01  is  set  to  FFH  to  indicate  the  start  of  the  control 
law  inputs  and  calculations.  The  port  is  set  to  00  to  mark  the  end 
of  the  calculations.  The  time  of  the  calculations  can  be  measured 
by  connecting  an  oscilloscope  to  this  output  port  and  observing 
the  pulse  length.  The  routine  INTCT  is  called  to  increment  a  coun¬ 
ter  that  keeps  track  of  interrupts.  The  routine  BLINK  is  called  to 
flash  the  instrument  panel  light  at  one-tenth  the  sampling  rate. 

There  is  basically  only  one  restriction  on  the  Flight 
Control  routine  for  the  proper  operation  of  CAS-1.  The  routine  must 
complete  execution  within  the  sampling  interval  or  no  other  routines 
will  run.  In  such  a  case,  there  would  be  no  available  time  for  key¬ 
board  inputs,  display  outputs,  or  memory  checks.  See  the  CAS-1 
source  listing  (not  included  in  the  thesis)  for  detailed  examples 
of  using  mode  set-up  and  flight  control  interrupt  service  routines. 
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